Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 12 April 2013 (MN KTgX style file v2.2) 



Star Formation & Dust Heating in the FIR Bright Sources of M83 



en : 

o 

Or 
< 

O 



O 

u 

Oh: 
6 

CO 

C3 



K. Foyle/f G. Natale, 2 f C. D. Wilson, 1 C. C. Popescu, 2 ' 3 M. Baes, 4 G. J. Bendo, 5 
M. Boquien, 6 A. Boselli, 6 A. Cooray, 7 D. Cormier, 8 I. De Looze, 4 J. Fischera, 9 O. L. 
Karczewski, 10 V. Lebouteiller, 8 S. Madden, 8 M. Pereira-Santaella, 11 M. W. L. Smith, 12 
L. Spinoglio, 13 R. J. Tuffs 3 

1 Dept. of Physics & Astronomy, McMaster University, Hamilton, Ontario, L8S4M1, Canada 
2 Jeremiah Horrocks Institute, University of Central Lancashire, Preston, PR1 2HE, UK 
3 Max Planck Institute fur KernPhysik, Saupfercheckweg 1, D-691 17 Heidelberg, Germany 
4 Sterrenkundig Observatorium, Universiteit Gent, Krijgslaan 281 S9, B-9000 Gent, Belgium 

5 UK ALMA Regional Centre Node, Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, University of Manchester, Oxford Road, Manch- 
ester Ml 3 9PL, United Kingdom 
6 Aix Marseille Universite, CNRS, LAM (Laboratoire d'Astrophysique de Marseille) UMR 7326, 13388, Marseille, France 

7 Department of Physics & Astronomy, University of California, Irvine,CA 92697, USA 

8 CEA, Laboratoire AIM, Irfu/SAp, Orme des Merisiers, F-91191 Gif-sur-Yvette, France 

9 'Canadian Institute for Theoretical Astrophysics, University of Toronto, 60 Saint George Street, Toronto, ON, M5S 3H8, Canada 
10 Department of Physics and Astronomy, University College London, Gower Street, London, WC1E6BT, UK 
11 Istituto di Astrofisica e Planetologia Spaziali, INAF, Via Fosso del Cavaliere 100, 1-00133 Roma, Italy 
12 School of Physics & Astronomy, Cardiff University, Queen Buildings, The Parade, Cardiff CF24 3A 
ls Istituto di Fisica dello Spazio Interplanetario, INAF, Via del Fosdso del Cavaliere 100, 1-00133 Roma, Italy 



accepted to MNRAS 



> 

oo 

(N 

m 

o 



% 



ABSTRACT 

We investigate star formation and dust heating in the compact FIR bright sources detected in 
the Herschel maps of M83. We use the source extraction code GETSOURCES to detect and 
extract sources in the FIR, as well as their photometry in the MIR and Ha. By performing in- 
frared SED fitting and applying an Ha based star formation rate (SFR) calibration, we derive 
the dust masses and temperatures, SFRs, gas masses and star formation efficiencies (SFEs). 
The detected sources lie exclusively on the spiral arms and represent giant molecular associ- 
ations (GMAs), with gas masses and sizes of 1O 6 -1O 8 M and 200-300 pc, respectively. The 
inferred parameters show little to no radial dependence and there is only a weak correlation 
between the SFRs and gas masses, which suggests that more massive clouds are less efficient 
at forming stars. Dust heating is mainly due to local star formation. However, although the 
sources are not optically thick, the total intrinsic young stellar population luminosity can al- 
most completely account for the dust luminosity. This suggests that other radiation sources 
contribute to the dust heating as well and approximately compensate for the unabsorbed frac- 
tion of UV light. 

Key words: galaxies: individual — , galaxies: ISM, galaxies: spiral, ISM: dust 



1 INTRODUCTION 

The Herschel Space Observatory (Pilbratt et al. 2010) has provided 
us with high sensitivity and angular resolution maps of nearby 
galaxies in the far-infrared, which has allowed us to spatially re- 
solve their cold dust emission. Dust plays a key role in the chem- 
istry of the interstellar medium, acting as a catalyst for the forma- 
tion of molecular gas, the fuel for star formation. However, it com- 
plicates our view of galaxies by obscuring UV and optical photons 
from stars and then re-radiating this light in the infrared. Thus, over 
a third of a galaxy's bolometric luminosity comes to us at these 
longer wavelengths (e.g. Draine et al. 2003; Bernstein et al. 2002). 



U In the past it has been very difficult to detect the dust emission in 
galaxies between 200 and 850/im. Ground-based telescopes have 
lacked the sensitivity and, prior to Herschel, space telescopes could 
not make detections at these wavelengths. 

Apart from the measurement of dust luminosities and masses, 
the maps from Herschel, spanning 70 to 500 fim, can be used to 
estimate two important quantities: the gas mass and the average in- 
tensity of the radiation field heating dust within galaxies. In the first 



f The first two lead authors have been co-equal contributors to the majority 
of the work presented in this paper. 
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case, since dust is generally well-mixed with gas, we can use dust 
mass estimates inferred by SED fitting as a proxy for the total gas 
mass, assuming a gas-to-dust ratio (e.g. Hildebrand 1983, Boselli 
et al. 2002, Eales et al. 2010, Eales et al. 2012). While atomic gas 
measurements are relatively well-known due to 21 cm line mea- 
surements, molecular gas measures are more challenging since we 
are forced to use an alternative tracer like CO rather than directly 
measure molecular hydrogen. This requires the calibration of this 
tracer, which is known to vary with environment, metallicity and 
density (e.g. Shetty et al. 2011). Herschel maps have been used 
to help spatially resolve this calibration factor (i.e. the X-factor) 
improving our view of the molecular gas component in galaxies 
(Sandstrom et al. sub. ). Typically dust emission maps also have 
higher resolutions, at least at the shorter wavelengths, and superior 
sensitivity to CO maps. 

Dust emission SED fitting also provides a measure of the aver- 
age radiation field energy density heating the dust or, alternatively, 
the average dust temperature. These latter quantities are connected 
to the luminosity of the heating sources and are therefore, in prin- 
ciple, useful to understand which radiation sources are heating the 
dust. In particular, in order to use dust emission to infer star for- 
mation rates, one would like to quantify the fraction of dust heating 
due to recent star formation. However, the dust can also be powered 
by an older stellar population (e.g. Popescu & Tuffs 2002, Groves 
et al. 2012, Bendo et al. 2012, Smith et al. 2012). 

Herschel has improved our view of the morphological struc- 
ture of the dust emission and has allowed us to better separate 
compact sources from diffuse emission in some of the far-infrared 
bands. We find both clumpy, compact regions associated with sites 
of recent star formation, detected also in Ha emission, and a 
smoother, more diffuse component (e.g. Verley et al. 2010). Al- 
though Herschel images have a resolution of 6" at the shortest 
wavelengths (70/im), at longer wavelengths it increases to 36" 
at 500 fim. Therefore, to date most studies have been forced to 
degrade the high-resolution PACS (Photodetector Array Camera; 
Poglitsch et al. 2010) maps (70 and 160 pirn) in order to match 
them to the longer wavelength SPIRE (Spectral Photometric Imag- 
ing REceiver; Griffin et al. 2010) data (i.e. Bendo et al. 2010, Smith 
et al. 2012, Bendo et al. 2012, Boquien et al. 2012, Aniano et al. 
2012). This has meant that these studies have averaged over large 
physical areas (i.e. pixel-by-pixel analysis at 36" resolution or az- 
imuthal averages). Thus, the structure evident in the PACS maps 
is mostly lost. In this work, we attempt to capitalize on the high- 
est resolution PACS maps, while still making use of the lower res- 
olution SPIRE maps, in order to study the compact, clumpy re- 
gions in the FIR continuum emission. Our method relies on the 
multi-wavelength extraction tool, GETSOURCES (Men'shchikov et 
al. 2012), which allows us to preserve the native resolution of the 
images. 

This paper uses this new technique on M83, a galaxy from the 
Very Nearby Galaxies Survey (PI: CD. Wilson). The close prox- 
imity of M83 (4.5 Mpc; Thim et al. 2003) affords us high spatial 
resolution (130 pc at 70 /im). M83 also has a strong spiral struc- 
ture and prominent bar, allowing us to investigate different regions 
within the galaxy. 

The main goal of this work is the investigation of the prop- 
erties of the compact sources detected in the FIR, in terms of gas 
masses, star formation and dust heating. Specifically, we address 
the following questions: 1) Which range of giant molecular cloud 
masses are associated with the FIR compact sources in M83?; 2) 
How efficient is the star formation within them?; and 3) Is the local 



star formation the predominant radiation source powering the dust 
emission? 

To this purpose, we developed a procedure consisting of the 
following steps. We use GETSOURCES to detect and extract com- 
pact regions from the FIR (70-350 pirn), which we call "clumps". 
We also measure the flux in the MIR (8 and 24 /im) and Ha using 
the clump area and position. In all cases a local background emis- 
sion component is subtracted. For each well-detected source, the 
corresponding dust luminosity, dust mass, average radiation field 
energy density and dust temperature are determined using a two- 
component dust emission SED fitting method of the MIR and FIR. 
Using the dust mass we infer a gas mass for each source, using a 
constant gas-to-dust ratio. We also derive the clump SFRs by ap- 
plying the Calzetti et al. (2007) calibration on the measured Ha 
and 24 /im fluxes. 

We structure the paper in the following sections: background 
and motivation (§2); observations and data reduction (§3); compact 
source extraction (§4); dust emission SED fitting (§5); star forma- 
tion rates and gas mass estimations; presentation of results (§7) and 
discussion (§8). We conclude with a summary of the main findings 
(§9). 



2 BACKGROUND & MOTIVATIONS 

2.1 Gas & Star Formation on sub-kpc scales 

During the last few years, high-resolution maps of nearby galax- 
ies in several wavelengths including UV, IR, emission line and ra- 
dio have become gradually available and opened the possibility to 
study the relation between gas mass surface density and star for- 
mation, known as Schmidt-Kennicutt law (hereafter S-K relation; 
Kennicutt 1998; Schmidt 1959) on scales of kpc/sub-kpc within 
galaxies. The S-K relation is commonly expressed as 



:AE£ 



(1) 



The value of the exponent N is particularly important and remains 
debated. N=l implies that the star formation efficiency (SFE) or gas 
depletion timescale is constant, while N> 1 implies that the SFE is 
greater for high-density regions, or that the timescale is shorter for 
high-density regions. While global averages across galaxies have 
shown that the total gas mass (HI + H2) surface density has a power 
law relation with the SFR surface density (i.e. N=1.4; Kennicutt et 
al. 1989), spatially resolved studies of star formation have shown 
that the SFR surface density is more tightly coupled to the molec- 
ular gas surface density (e.g. Wong & Blitz 2002 and Bigiel et al. 
2008). 

Reliable measurements of SFR within galaxies can be per- 
formed by combining optical/UV and infrared flux measurements, 
according to calibrations such as Calzetti et al. (2007) or Zhu et 
al. (2008) for Ha+24/im and Bigiel et al. (2008) for UV+24/^m. 
In these empirically based calibrations, dust emission is used to 
estimate the fraction of optical/UV light attenuated by dust. These 
calibrations are now considered the most accurate way to determine 
SFRs within galaxies, at least when used on galaxy regions which 
are bright enough (SFR >= 0.001 M© yr _1 , Kennicutt & Evans 
2012). 

Spatially resolved studies of the SFR typically involve aper- 
ture photometry of the SFR tracers. The size of the aperture has 
been primarily dictated by the beam sizes. Typically background 
emission is subtracted from the SFR tracers either locally, in an- 
nuli directly surrounding the apertures, or, in the case of crowded 
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regions, more extended zones are used (Calzetti et al. 2005). Back- 
ground subtraction is typically performed because the star forma- 
tion tracers have contamination from other sources not directly as- 
sociated with star formation 1 !! The 24/im emission includes cirrus 
emission from an older stellar population (Popescu & Tuffs 2002), 
which may contribute as much as 30% of the emission (Kennicutt 
et al. 2007). Ha emission may also include diffuse ionized gas (Fer- 
guson et al. 1996). Background subtraction has proven to have ef- 
fects on the S-K relation. Subtracting a diffuse component more 
readily suppresses the fainter star forming regions (Liu et al. 201 1). 
Thus, studies that have used aperture photometry on star forming 
regions and background subtraction have found a super-linear rela- 
tion (e.g. Kennicutt et al. 2007). Meanwhile pixel-by-pixel analyses 
that have not subtracted a diffuse component have recovered a lin- 
ear relation (e.g. Bigiel et al. 2008). 

The measures of star formation tracers are quite advanced due 
to the relatively high resolutions available but tracing gas in extra- 
galactic studies has proven to be more complicated. While mea- 
surements of the atomic gas component are available by using the 
21 cm line emission (e.g. THINGS; Walter et al. 2008), cold molec- 
ular hydrogen has proven more challenging. Due to the low mass of 
the molecule, it requires high temperatures to excite the rotational 
transitions of molecular hydrogen. Thus, it is virtually impossible 
to measure the total amount of molecular hydrogen directly. In- 
stead, an alternative tracer like CO is commonly used (Bolatto et al. 
2013). Typically conversions from CO are done using a constant X- 
factor. However, this factor can vary considerably depending upon 
the metallicity, density and temperature (Kennicutt & Evans 2012). 
Accurate measures of the molecular component associated with 
star formation are particularly important, because it seems that star 
formation is more directly coupled to molecular gas than to total 
gas. 

Due to the comparatively low resolution and sensitivity of the 
molecular gas maps, only few studies (e.g. Rahman et al. 2011) 
have subtracted diffuse emission from the gas. However, diffuse gas 
unrelated to star formation sites is most likely present. Low mass 
clouds that are unresolved will present themselves as a diffuse com- 
ponent and may not host star formation. Ideally, one should also ac- 
count for a diffuse gas component. Furthermore, there is growing 
evidence that low mass GMCs may not have a corresponding star 
forming region (e.g. Hirota et al. 201 1). 

In light of the issues surrounding the molecular gas measures 
we are motivated to capitalize on the high resolution dust emission 
maps from Herschel. By employing a gas-to-dust ratio we present 
an alternative way to probe the star formation efficiency on spa- 
tially resolved scales. Inspired by the aperture photometry studies 
of star forming regions, we employ the dust maps to isolate and ex- 
tract clumpy regions which are likely associated with individual or 
groups of molecular clouds. For both the star formation tracers and 
the dust we can also subtract a diffuse background component, thus 
treating both the star formation and gas tracers in a similar fashion. 



2.2 Dust heating within galaxies 

In general, dust emission within galaxies without AGNs is powered 
by radiation coming from both sites of recent star formation and 



x The relative amount of background emission that needs to be removed is 
still the subject of much debate. See e.g. Leroy et al. (2012) for a discussion 
on the nature of the diffuse MIR emission. 



from more evolved stellar populations. However, there is a long- 
standing debate about the exact fraction of dust heating contributed 
by each stellar population (e.g. Law et al. 201 1, Boquien et al. 201 1, 
Bendo et al. 2012, for recent references), which depends on several 
factors: the intrinsic emission spectral energy distribution of the 
stellar populations, the dust mass and optical properties, the rela- 
tive dust-star geometry. 

Recent observational works investigating the origin of the ra- 
diation heating the dust in nearby galaxies have looked for corre- 
lations between the source dust temperature (or, alternatively, FIR 
colour) and source stellar populations luminosities, as traced by 
SFR for the young stars and NIR luminosity for the old stars (e.g. 
Bendo et al. 2012, Boquien et al. 2011, Smith et al. 2012, Foyle et 
al. 2O12). 2 

The presence of an observed correlation of this kind has of- 
ten been used as evidence that a particular stellar population is the 
dominant source heating the dust. Although it might seem intuitive 
that a higher dust temperature should correspond to a stellar popu- 
lation with a higher intrinsic luminosity, deducing the origin of the 
radiation heating the dust in this way can be potentially misleading. 
The reason is that the intrinsic luminosity of a stellar population 
and the average dust temperature are not necessarily proportional 
to each other, even in the case where the stellar population in ques- 
tion is the dominant source of dust heating. 

The luminosity of a stellar population gives the amount of to- 
tal radiative energy per unit time injected by a stellar population lo- 
cally in the interstellar medium. Therefore, one can expect at most 
a proportionality between the observed dust luminosity and the in- 
trinsic stellar population luminosity for a set of sources, provided 
all the sources are mainly heated locally by the same kind of stellar 
population and heating from radiation sources external to the areas 
considered is negligible^. 

Specifically, one can express the total dust luminosity Ld us t 
for a source as: 



-^dust — -^stars (I 



(2) 



where Lstars is the total intrinsic luminosity from a certain stellar 
population and r is the source luminosity-weighted optical depth. 
The optical depth can be expressed as: 



r out 

-^stars 
Tint 



(3) 



where L°^ IS is the escaping unabsorbed source stellar population 
luminosity. All the aforementioned luminosities are integrated over 
wavelength spanning the entire emission spectral range. From Eq. 
[2]one can see that, if r has similar values for a sample of sources, a 
correlation between Ldust and Lstars will be found. Note also that 
Eq[2] implies that I/dust cannot exceed L"t* rs . If this is observed 
for a sample of sources, it would mean that the dust is significantly 
heated by radiation coming from another stellar population beyond 
that associated with L"ta rs ■ 

In order to have a proportionality between the source intrin- 
sic stellar luminosity and the source average dust temperature, an 

2 The quantities typically considered are surface densities of SFR and 
NIR luminosity rather than total values. For pixel-by-pixel analyses this 
amounts to multiplication with a constant, thus it does not affect the nature 
of the correlation with dust temperature. 

3 The latter assumption can easily break down when considering arbitrary 
galactic regions. For the compact FIR sources associated with star formation 
regions, the heating is usually thought to be dominated by the local young 
stellar populations but see further discussion in §8.3. 
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additional assumption is required: each of the detected sources 
should be associated with a similar amount of dust mass, Afdust ■ 
Assuming that the dust emission can be described by a modi- 
fied blackbody function such that k„.B„ (Tdust), where k„ oc v^ 
is the absorption coefficient, /3 is the emissivity, B v the Planck 
function and Td ust the dust temperature, it can be shown that 
Ldust oc Mdust / KvB v (Tdust)dv oc MdustT^+f. By combining 
the latter relation with Eq.|2] it follows that: 



T int 

-^stars 

A/d„ Rt , 



T 4+/3 
J dust 

(1-e-) 



(4) 



If Mdust and r have similar values for all the sources con- 
sidered, the average intensity of the radiation heating the dust ap- 
proximately scales only with the stellar luminosity and, as a conse- 
quence, the dust temperature Td us t is directly related only to L^ IS 
through Eq. [4] In this case, it is likely that a specific stellar popu- 
lation is heating the dust, if one finds that the intrinsic luminosity 
of that stellar population correlates with the observed dust temper- 
ature. However, the amount of dust mass, the dust-star geometry 
and, therefore, the value of r, can be substantially different for each 
source. Thus, in general, it cannot be expected that the luminosity 
of a stellar population is correlated with the dust temperature, even 
if that stellar population is responsible for the dust heating. 

Eq. [2] and [4] are almost equivalent to each other and they both 
express the relation between the dust luminosity and stellar popula- 
tion luminosity. The only difference is that Ldust in Eq.f2]includes 
all the dust and PAH emission throughout the entire infrared range. 
In Eq. [4] we assumed that the dust emission can be modelled by a 
single modified blackbody curve. This is a good approximation if 
one considers only the FIR region of the emission spectra. How- 
ever, usually the FIR luminosity is the dominant contribution to 

-^dust- 

Thus, Eq.|4]can be used to probe the extent to which a given 
stellar population is powering the dust emission. We can do so by 
comparing the local, cold dust temperature with a measure of the 
local stellar population luminosity divided by the corresponding 
dust mass. As before, if the value of r is varying in a relatively 
small range, a correlation between Lstars /A/dust and Td us t will 
be observed if the stellar population considered is responsible for 
heating the dust. 

For example, one can use the SFR value as a tracer of the UV 
luminosity of the young stellar populations associated with each 
source. Therefore, if the radiation from young stars dominates the 
dust heating, the SFR/Md us t ratio is expected to be more tightly 
coupled to the dust temperature than the SFR alone. Thus, an ob- 
served correlation between SFR/Afdust and Td us t would suggest 
that star formation is powering the observed dust emission. 

Given this, we are motivated to include a comparison of the 
measured SFR/ A/dust (or, equivalently, SFE if one assumes a con- 
stant dust-to-gas mass ratio) with Td us t for the sources we detect, 
in addition to simply a comparison of the SFR and dust tempera- 
ture, which is typically seen in the literature. The results of these 
comparisons are shown in Sect. §7 and discussed in Sect. §8. 



3 OBSERVATIONS 

In this work we use the Herschel far-infrared maps of M83 from 
the Very Nearby Galaxies Survey (PI: C. D. Wilson) as well as 
ancillary mid-infrared and Ha maps to trace dust/PAH emission 
and star formation rate. Fig.[T]shows M83 at each of the wavebands 
considered. 



3.1 Far-Infrared Images 

We use far-infrared (FIR) images from the Herschel Space Obser- 
vatory to trace cold dust emission. We use 70 and 160 pm maps 
taken with the Photodetector Array Camera (PACS; Poglitsch et al. 
2010) and 250 and 350 /jm maps taken with the Spectral Photo- 
metric Imaging REceiver (SPIRE; Griffin et al. 2010)f The PACS 
images are processed using both HIPE v5 and SCANAMORPHOS v8 
(Roussel 2012) and the SPIRE images are processed with HIPE and 
BRIGADE (Smith et al. 2012). The PACS images were corrected 
from the v5 photometric calibration files to v6 with corrective fac- 
tors of 1.119 and 1.174 for the 70- and 160- pm maps, respec- 
tively. The SPIRE images are multiplied by 0.9828 and 0.9839 for 
the 250 and 350 pm maps respectively in order to convert from 
monochromatic intensities of point sources to monochromatic ex- 
tended sources. The images are kept in their native resolution with a 
FWHMof thePSFof 6.0", 12.0", 18.2" and 24.5" for the 70, 160, 
250 and 350 pm maps, respectively. For more details on how the 
images were processed we refer to Bendo et al. (2012) and Foyle et 
al. (2012, hereafter F12). 



3.2 Mid-Infrared Images 

We trace the warm dust and PAH emission using mid-infrared maps 
(MIR) taken from the Spitzer Local Volume Legacy Survey (Dale 
et al. 2009). Specifically, we use the 8 and 24 pm maps from IRAC 
and MIPS instruments. We subtract the stellar component of the 
emission in the 8 pm map using a scaling of the IRAC 3.6 /jm map, 
according to the relation provided by Helou et al. (2004): F„(8 pm, 
dust) = -F„(8 ^m)-0.232_F„(3.6 pm). In the 24pm map the nuclear 
region is saturated, so this region is excluded from our analysis. 
Because of reasons explained in §4, we degrade the resolution of 
the MIR maps to 6" in order to match the resolution of the 70 ^m 
map. 



3.3 Ha Images 

We use continuum subtracted Ha maps from the Survey for Ioniza- 
tion in Neutral Gas Galaxies (SINGG; Meurer et al. 2006). We cor- 
rect the Ha maps for Galactic extinction using a factor 1.167 from 
the NASA/IPAC Extragalactic Database that is based on Schlegel 
et al. (1998). As for the MIR maps, we degrade the Ha map to a 
resolution of 6" matching that of the 70 pm map. We use the Ha 
maps in conjunction with the 24 pm MIR map, in order to mea- 
sure the star formation rate (SFR) of the extracted compact sources 
(discussed in greater detail in §6). 



4 COMPACT SOURCE EXTRACTION 

In order to extract compact sources in the FIR maps (70-350 pm), 
we use the new multi-scale, multi-wavelength tool, GETSOURCES 
(Men'shchikov et al. 2012). GETSOURCES is specifically designed 
to work with the FIR images of Herschel. The data from Herschel 
spans a range of angular resolutions from 6.0" to 36.0" and, thus, 



4 We performed an initial test including the 500 pm map as well, but we 
found that the uncertainties on the source fluxes were so large that they 
were not useful to constrain the source SEDs. This was caused by the low 
resolution of the 500 pm map (PSF FWHM=36"). Thus, we decided to not 
include the 500 pm map in our analysis. 
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Figure 1. M83 at different wavelengths at the resolution used in this study denoted by the white circle in the bottom right of each panel. Upper row: Ha from 
SINGG survey and Spitzer IRAC 8 and MIPS 24 fiva maps. Lower row: Herschel 70, 160, 250 and 350 /im maps. 



any source extraction code must be able to handle these extremes, 
which poses concerns for source blending at longer wavelengths. 

Rather than extracting sources directly from observed images 
(i.e. GAUSSCLUMP; Stutzki & Gusten 1990, CLUMPFIND; Williams 
et al. 1994 or SEXTRACTOR; Bertin & Arnouts 1996), GET- 
SOURCES analyzes spatial decompositions of the images across dif- 
ferent scales and different wavelengths. Wavelength independent 
images are generated to detect sources at each spatial scale. The 
original images are then used to perform photometry on the de- 
tected sources. This procedure takes into account the different an- 
gular resolutions of the maps, background subtraction and source 
blending. We briefly describe the key steps of the process here, but 
refer the reader to Men'shchikov et al. (2012) for a detailed descrip- 
tion on how GETSOURCES extracts and measures the properties of 
compact sources, including a comparison with other similar codes. 

After the images are aligned to the same spatial grid, GET- 
SOURCES decomposes of the original maps into single scale detec- 
tion images. This is done by using a process of successive unsharp 
masking, where the original images are convolved with Gaussians 
and subtracted successively, in order to enhance the visibility of 
emission on different scales. The FWHM of the Gaussians varies 
between twice the pixel size to a maximum of 18 times the resolu- 
tion of the image or the image size. The image resolution is the only 
information the user needs to provide. The background and noise 
in the single scale detection images is then removed by intensity 
thresholding. The clean single-scale detection images at each wave- 
length are then combined into single-scale wavelength-independent 
detection images, allowing for the use of all the information across 
all bands simultaneously for the source detection. Many detection 
codes rely on independent catalogues at each waveband which are 
then matched using an association radius, which can introduce large 
unknown errors. By dividing the images into single scales, this pro- 
cess can be avoided. 

On the combined single-scale detection images, a given source 



will appear at a small scale and gradually get brighter until it is seen 
at a scale roughly the true size of the source. Beyond this, the source 
begins to vanish again. GETSOURCES tracks the evolution of the 
source through the spatial scales and creates source masks to iden- 
tify the sources. The scale where the source is brightest provides an 
initial estimate of the source footprint. A source must have a signal- 
to-noise ratio of at least 3a in at least two bands in order to be con- 
sidered detected. Coordinates of the sources are determined using 
the moments of the intensities. Once the sources have been detected 
in the combined images over all spatial scales, GETSOURCES per- 
forms source flux measurements on the observed images at each 
waveband and simultaneously subtracts a background by interpo- 
lating under the sources footprints. Partially overlapping sources 
are "deblended" using an iterative process. 



Upon completing the extraction and measurements, the user is 
supplied with a table giving the location and size of each source 
at each waveband. Due to the fact that the resolution decreases 
with increasing wavelength, the apparent source size increases with 
wavelength. The measured properties also include total flux, peak 
flux, degree of source blending, monochromatic and global detec- 
tion significance. The detection significance at each wavelength is 
essentially a signal-to-noise ratio, which is determined by the ra- 
tio of the peak flux of the source to the standard deviation in an 
annulus surrounding the source on the detection maps. The global 
detection significance is determined by the square root of the sum 
of the squares of the detection significance at each wavelength. As 
discussed in the following sections, we apply a minimum threshold 
for the global detection significance value to remove some sources 
not well detected. 



We now briefly describe the parameters we defined in the 
source extraction process. 
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4.1 Preparing Images 

We supply all images (Ha, MIR and FIR) to GETSOURCES, but 
only use the FIR images to detect the compact sources, since we 
seek to trace the location of cold compact clouds. 

The images are aligned to the same grid and converted to 
MIy sr" 1 , with a pixel size of 1.4". For the SPIRE images, GET- 
SOURCES uses the beam areas in order to convert the units from 
Jy/beam (423 ± 3, 751 ± 4, 1587 ± 9 arcsec 2 , for the 250, 350 
and 500 /im images respectively). The images are all aligned to the 
WCS of the 70/im maps. Observational masks are created, which 
denote the image area over which GETSOURCES is meant to look 
for sources, which speeds the detection process. The alignment and 
masks are visually checked. 

Since GETSOURCES has been designed to work with FIR Her- 
schel maps, we decided to degrade the MIR and Ho maps to the 
angular resolution of the 70 ^m map. This should reduce possible 
systematic effects due to the use of the program on a larger range 
of spatial resolutions than the one on which it has been designed 
and tested. Furthermore, we performed a test to examine the ef- 
ficiency of GETSOURCES in recovering source fluxes at different 
resolutions, by running the program on a set of convolved 70/im 
maps. In this test, source fluxes are expected to be the same at each 
resolution. We found the accuracy of the source flux measurement 
decreases with poorer resolution, but the average systematic effect 
does not seem to be large enough to affect our results substantially 
(see Appendix A for more details). 



4.2 Compact Source Measurements 

The source detection has been performed using only the FIR maps, 
since we aimed to select sources due to their FIR brightness and 
not necessarily MIR and Ha counterparts. However, MIR and Ha 
maps have been provided to GETSOURCES for source photometry. 
GETSOURCES detects 186 compact sources across M83. However, 
in the analysis described in the following sections, we only include 
those sources with a detection significance above 20. Furthermore, 
since a high global detection significance does not necessarily im- 
ply that the photometry is accurate for each single waveband, we 
further removed all the sources with a photometric signal-to-noise 
ratio lower than Ict in at least one band between 8 and 250 /im (that 
is, we occasionally retained sources detected in all bands but the 
350 fim band). After this source selection, our sample is reduced to 
121 sources. 

Fig. [2] shows the location of all the detected sources overlaid 
on the 70 /im map, with colours denoting the detection significance. 
Sources that are retained in the study are denoted with triangles, 
while sources that are excluded are denoted with squares. We note 
that the sources lie almost exclusively on the bar ends and the spiral 
arms. Very few sources are detected in the interarm regions and 
along the bar. This does not mean that compact sources do not exist 
in these regions, but rather, the sources there are not bright enough 
to be detected. Thus, it seems that only the bar ends and spiral arms 
harbour bright compact sources that we can detect. 

Our spatial resolution at 70 /im allows us to measure intrinsic 
sources sizes only for sources having radii larger than 130 pc. We 
find a median radius of 150 pc at this waveband. This implies that 
the majority of our sources are unresolved or barely resolved. Thus, 
it is not possible to determine a physical radius for all the sources. 

At each waveband we obtain the total flux of each of the 
extracted sources (hereafter, F SOU rce, hb) and of the corresponding 
background emission, which has been interpolated and subtracted 
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Figure 3. Median value of the ratio of the measured source flux to the total 
flux in the source footprint for both the high (full circles) and low (open 
squares) source background subtractions. At longer wavelengths the back- 
ground is greater due to the increasing size of the footprint. 



by GETSOURCES. Because GETSOURCES interpolates the back- 
ground from regions very close to the sources, this background 
measurement is actually determined in the vicinity of the spiral 
arms. We find that the background makes up more than 50% of the 
flux in the footprint. Fig. [3] shows the relative flux of the source 
to the total flux in the footprint. The relative fraction decreases 
with wavelength mainly because, due to the poorer resolution at 
longer wavelengths, the footprint size increases and thus a rela- 
tively larger amount of background flux is included in the footprint. 
However, part of this relative flux variation might be also due to an 
intrinsic higher fraction of diffuse background emission at longer 
wavelengths. We find that at all wavelengths the source makes up 
less than 50% of the emission in the footprint. For this reason, we 
refer to this way of performing the source photometry as "high- 
background" source flux measurement. In order to identify possi- 
ble systematic uncertainties due to the way the background is sub- 
tracted, we also consider an alternative measure for the background, 
estimated in the interarm regions, as shown in the next subsection. 



4.3 Low-Background Source Flux Measurement 

The background estimated by GETSOURCES is basically a measure 
of the smooth emission component associated mainly with the spi- 
ral arms. However, typically the photometry of star forming regions 
is done by subtracting an average brightness measured in the inter- 
arm regions close to the sources (i.e. Calzetti et al. 2007 and Liu 
et al. 201 1) and, as said before, it is important to check if there are 
substantial differences on the final results if the background mea- 
surement is performed differently. 

For these reasons, we decided to perform aperture photometry 
in several interarm regions selected by eye (see Fig. |2]- Specifi- 
cally, we selected 18 interarm regions, which are all sufficiently far 
away from the detected sources at all wavelengths. The diameter 
size of the apertures is 25" which is large enough to obtain a good 
estimate of the rather uniform interarm emission. For each aperture 
we measured the average surface brightness with an uncertainty de- 
rived from the standard deviation of the surface brightness within 
each aperture and the flux calibration uncertainty. 

Once we obtained the flux measurements for the interarm 
apertures, we performed a "low background" source flux measure- 
ment in the following way: for each source, we consider the flux in 
the full footprint determined by GETSOURCES and subtract the in- 
terarm surface brightness, estimated from the aperture that is clos- 
est to the source, multiplied by the area of the source footprint (typ- 
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Figure 2. M83 at 70 fim with the compact source locations marked with colors corresponding to the detection significance (Log 10 scale). A detection 
significance of 20 is used as a cut-off for retaining sources for further analysis. Sources that are used in the study are denoted with triangles and excluded 
sources are denoted with squares (see text for criteria). Apertures showing the measurement locations of interarm emission are displayed in white. 



ically, the distance between the centres of each source and the clos- 
est interarm region is about 1 kpc or less)^. Since the GETSOURCES 
background is usually larger than that determined in the interarm 
region, in general the "low background" (LB) source fluxes are 
higher than the "high-background" (HB) ones. Fig.[3]shows that the 
LB measurements have higher fluxes relative to their background 
at all wavelengths than HB flux measurements. 



5 SED FITTING 

In the absence of a strong contribution from AGNs, the galaxy 
dust/PAH emission SED on sub-kpc/kpc scales can be modelled as 
the sum of two emission components: a warm component largely 
emitting in the MIR, produced by dust in photo-dissociation re- 
gions (PDRs) and heated predominantly by young stellar popula- 
tions; and a diffuse emitting component, emitting mostly in the 
FIR and MIR PAH line emission, which can be powered both by 
older stellar populations and by the fraction of UV photons es- 
caping from PDRs. This concept is at the base of the SED fitting 
method developed by Natale et al. (2010, NA10), which we used in 
this work to fit all the well-detected source emission SEDs derived 
in §4. In the following we explain the main features of this SED 
fitting method, together with minor updates, and its application to 
our data set. We refer the reader to NA10 for additional details. 



5 We note that there are variety of possible approaches for subtracting 
background emission. Here we contrast two techniques, one which utilizes 
an interpolation scheme in the vicinity of the source and one that relies on 
a more distant estimate in the interarm regions. 



Using the mentioned fitting method, each observed source 
dust emission SED is fit by combining two infrared SED com- 
ponents (see Fig|4j. The first component is a PDR SED template 
which has been selected by Popescu et al. (2011) among the mod- 
els of Groves et al. (2008) because it provides a good fit to the dust 
emission from Milky Way star formation regions (specifically the 
chosen model is for compactness parameter log(C) — 6.5, solar 
metallicity and hydrogen column density log(N) = 22, see §2.8 
of Popescu et al. 2011 for more details). The second component, 
suitable to fit the diffuse dust emission, is taken from a grid of uni- 
formly heated dust emission templates, obtained by using the dust 
emission code of Fischera & Dopita (2008). The diffuse dust emis- 
sion is calculated assuming a Milky Way dust+PAH composition 
(exact dust model parameters can be found in Table 2 of Fischera 
& Dopita 2008) and taking into account the stochastic heating of 
grains following the method of Guhathakurta & Draine (1989), 
combined with the step wise analytical solution of Voit (1991). 
Since we assumed a fixed dust/PAH composition, each element of 
the grid is determined only by the parameters of the radiation field 
heating the dust. The adopted spectral shape of the radiation field is 
the classical Mathis et al. (1983) profile, which was derived for the 
local interstellar radiation field, scaled by two linear factors: Xuv, 
which multiplies the whole curve, and Xcoi which multiplies only 
the optical part of the Mathis spectra (see appendix B.2 of NA10). 
Therefore, xuv can be seen as the intensity of the UV radiation 
field and Xcoi as the optical to UV ratio in the units of the Mathis 
etal. (1983) profile. 

Compared to NA10, we extended the size of the grid of the 
diffuse dust templates, in order to cover a larger range of possi- 
ble radiation field parameters. Specifically, xuv and Xcoi are both 
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allowed to vary between 0.1 and 10, a range which is reasonably 
large to include all the plausible values of the diffuse radiation field 
intensity and colour within galaxies (note that even higher radia- 
tion fields and, therefore, warmer dust are associated with the PDR 
component in our SED fitting procedure). We also point out that 
there is a degeneracy between Xcoi and the relative abundance of 
PAH and solid dust grains, in the sense that they both affect the 
8 /mi/FIR ratio (see §B. 2 of NA10). 

One of the main differences between our dust emission models 
and those of Draine & Li (2007) is that in their models the PAH 
abundance is varied but the optical/UV intensity ratio is fixed to 
the Mathis et al. (1983) value. Therefore, this should be taken into 
account when considering results involving the Xcoi parameter. 

We performed a x 2 minimization fitting of the data with the 
two dust emission components by varying four free parameters: 
Xuv and Xcoi, the radiation field parameters defined above; Md 
the dust mass of the diffuse dust component; F24, the fraction of 
24 /im emission provided by the PDR component. For each pair of 
Xuv and Xcoi values, the parameters F24 and Md us t correspond to 
linear scaling factors for the amplitudes of the PDR template and 
the diffuse dust component respectively. The fit takes into account 
colour corrections calculated for each SED template following the 
conventions adopted for each instrument (See also Sect. 5 of NA10 
for details and the Spitzer and Herschel observer manuals for the 
colour correction definitions). 

The uncertainties on each parameter are calculated by analyz- 
ing the variation of \ 2 around the minimum found by the fitting 
procedure. The one sigma interval is defined as the minimum vari- 
ation of a given fitting parameter around its best fit value, which 
produces a variation A\ 2 = X 2 ~Xmin always higher than 1, inde- 
pendently of all the possible values of the other fitting parameters. 
The program checks also if "islands" of \ 2 l° w values (such that 
Ax 2 < 1) are present, which are detached from the region where 
the minimum of x 2 has been found. In that case, a conservative un- 
certainty covering all the regions of low x 2 values is provided for 
the fitting parameter. 

Apart from the fitting parameters, our SED fitting procedure 
provides the total dust emission luminosity and the luminosities of 
each SED component. However, in contrast to blackbody fits, our 
method does not provide a single average dust temperature, since 
the output total spectra of the dust emission code is determined by 
probability distributions of dust temperatures, which are different 
for each dust grain size and composition. Furthermore, in princi- 
ple different combinations of UV and optical radiation field energy 
densities can cause the dust to have approximately the same aver- 
age cold dust temperature (that is, similar FIR peak wavelength), 
which would not be immediately evident by comparing different 
pairs of xuv and Xcoi values. However, for the purposes of com- 
parisons with other works that have relied on blackbody fits, it is 
useful to quote an average dust temperature. Taking advantage of 
the fact that the cold FIR part of the diffuse dust emission compo- 
nent (A >70 /im) can be well fit with a modified blackbody curve 
with a dust emissivity index of /3=2 (consistent with observational 
results, i.e. F12, Davies et al. 2012, Auld et al. 2013), it is straight- 
forward to associate to each diffuse dust template the dust tempera- 
ture of the modified blackbody which best reproduces its FIR part. 
We will refer to this cold dust temperature as Tdust- In Appendix 
B we compare the dust masses and temperatures found using the 
SED fitting method outlined here to those determined directly from 
a modified blackbody function fit to the FIR wavelengths. 



6 STAR FORMATION RATES AND GAS MASSES 

Although the sources are detected in the FIR bands, measurements 
of the flux in the source footprint in the Ha and MIR maps are also 
made. Ha and 24/im emission can be used in conjunction to trace 
recent star formation - both unobscured and obscured (Calzetti et al. 
2007). GETSOURCES also interpolates and subtracts a background 
for these images as well. A visual comparison of the source mea- 
surements in Ha confirmed that the source footprints are centred 
on bright emission peaks. In this way, we feel confident that the 
flux measured in the footprints is directly related to an active star 
forming region. 

For each source, we combined the fluxes measured in Ha and 
24 /im emission to derive a star formation rate by using the calibra- 
tion of Calzetti et al. (2007): 

SFR[M Q yr _1 ] = 5.3 x 10~ 42 (L(H a ) + 0.031L(24fim)). (5) 

There are two key differences between our SFR and that de- 
rived by Calzetti et al. (2007). We detect the compact sources in the 
FIR maps rather than the Ha. This means we essentially estimate 
the SFR in the region of a FIR compact source. However, a vi- 
sual check shows that almost all sources in the FIR are also present 
in Ha. A second difference is related to how the background is 
treated. Here, we rely on the background subtraction performed 
by GETSOURCES. This background is interpolated in the region 
surrounding the sources and takes into consideration neighbouring 
sources and deblending. In Calzetti et al. (2007), the background 
is determined in 12 rectangular regions surrounding sources. These 
regions cover a large number of pixels and extend well beyond the 
local neighbourhood of the sources. The mode of these regions is 
then used as a measure of the background. Our method, produces 
a more 'local' background, which is greater than that which would 
be found by using the same method of Calzetti et al. (2007). In or- 
der to check the effect of choosing a different background level, 
we also determine source star formation rates by using the Ha and 
24 /im emission in the total source footprint, after subtraction of 
the background estimated from the nearest interarm aperture. This 
is exactly the same approach as for the LB measure for the compact 
source dust emission. 

We note a potential problem with the SFR as outlined above. 
During the source extraction, it is not known which sizes the 
sources will potentially have. The size of the region plays an im- 
portant role in deciding whether the calibration described above 
can be used. Regions which are too small may have SFRs which 
are too low. It is known that SFRs below ~ 0.001 Mq yr _1 can 
be problematic because there may be incomplete sampling of the 
IMF and the assumption of continuous star formation in the last 
few Myr may not be valid (Kennicutt & Evans 2012; Leroy et al. 
2012). Furthermore, the adopted calibration assumes that all the 
ionizing photons are absorbed by the gas in the HII regions. How- 
ever, a fraction of ionizing photons could escape the star formation 
regions before ionizing the gas or, instead, could be absorbed by 
dust (see e.g. Boselli et al. 2009, Calzetti 2012, Relano et al. 2012). 
Due to all these concerns, it is useful to consider the SFR measured 
in this process as a "corrected" Ha luminosity. Whether this lumi- 
nosity can accurately trace the SFR will depend on the source size 
and the value of the SFR. 

We estimate the source gas masses Af gas by making use of the 
fact that the dust mass can be used as a proxy for the total gas mass 
(both molecular and atomic), provided one can translate the dust 
mass with a constant gas-to-dust ratio (GDR) to a gas mass. Dust 
mass is known to be better correlated with the total gas mass than 
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Figure 4. Example of source flux measurement and SED fitting. The maps on the left show a small area of the disk of M83 around one of the detected sources. 
The ellipses overplotted on each map delineate the source footprint determined by GETSOURCES. On the right, we show the two-component SED fit of the 
same source with best fit parameter and SFR values listed aside. The dashed and the dotted line represent the diffuse dust emission component and the PDR 
emission component respectively. 



with the atomic or molecular component alone (i.e. Corbelli et al. 
2012). 

We assume a constant GDR of 100 for the compact regions. 
While, there is some evidence that the GDR may vary within galax- 
ies (e. g. Sandstrom et al. 2013), particularly with metallicity, given 
that our sources lie almost exclusively along the spiral arms, it is not 
obvious that a varying GDR should be used for these sources. For 
example, the metallicity values for M83 are based on azimuthal av- 
erages, which include interarm regions. Thus, employing a varying 
GDR based on metallicity values would require introducing another 
assumption, namely that the GDR depends more on radius than 
arm and interarm regions. In the pixel-by-pixel analysis of F12, 
the GDR was found to be relatively constant on the spiral arms (~ 
100) and, in contrast to the interarm regions, declined only slightly 
with radius with the inner regions of the spiral arm having a GDR 
of 130 and the outer tips of the spiral arm having a GDR of 120. 
Given the current uncertainties and possible systematic effects, we 
simply adopt a value of 100, which is in the range found by F12 
(mean value of 84 and standard deviation of 40). 

By combining the SFR and gas mass measurement for 
each source, we estimate the star formation efficiency defined as 
SFE=SFR/A/ g(ls for both types of compact source flux measure- 
ment. We stress that this SFE is different from that presented in 
other works (Leroy et al. 2008; Bigiel et al. 2008 etc.) which de- 
scribe the SFE using only the molecular gas component. However, 
it is likely that the dust is mostly tracing molecular gas in these 
regions. M83 is known to be a galaxy with a dominant molecular 
component and, over the region we consider, the molecular gas con- 
stitutes 80% of the total gas (Crosthwaite et al. 2002). F12 found 
that on the spiral arms, where most sources lie, the molecular gas 
component has column densities on the order of 100 Mq pc~ 2 , 
whereas the atomic gas had column densities of less than 10 Mq 
pc- 2 . 



7 RESULTS 



Upon completing the SED fitting, as described in §5, we further 
remove sources that have high (total) \ 2 values, that is % 2 > 10, 
or have a dust mass uncertainty greater than a factor of 2. In the 
first case we remove all the sources which are not well-fit by our 
model (mainly because the observed SED shape is too irregular to 



be fitted with our two-component model), while in the second we 
remove those sources with flux uncertainties that are too large at 
some wavebands and, therefore, cause the fitting parameters to not 
be constrained within small ranges. After applying this selection 
criteria, we are left with 90 sources. Appendix D presents tables of 
the extracted flux and SED fitting parameters for both the HB and 
LB source measurements. In this section we present the results for 
this set of sources. Specifically we describe the distributions of the 
inferred source parameters in§7.1, their radial variation in §7.2 and 
their interdependence in §7.3. 



7.1 Inferred Source Parameter distributions 

In order to show the distributions of inferred parameters and com- 
pare them for the two types of source flux measurements, we cre- 
ated box-and-whisker plots which are shown in Fig. [5] [6] and [7] 
The solid line in the boxes shows the median value and the box de- 
lineates the 25th and 75th quartiles. The lines extend to either the 
maximum and minimum values or to 1.5 times the 75th and 25th 
quartiles. If there are values beyond the later range, they are de- 
noted with open circles. In the following we describe in detail the 
distributions for each set of parameters, that is, the SED fitting pa- 
rameters, the source luminosities and the star formation rates and 
efficiencies. 



7.1.1 SED fitting parameters 

As explained in §5, the SED fitting has four free parameters Xuv, 
Xcoi, Mdust and F24,. The box-and-whisker plots of these parame- 
ters and of the dust temperature are shown in Fig|5] 

The xuv values (see bottom left panel), representing the UV 
radiation energy density in the units of the Mathis et al. (1983) pro- 
file (hereafter, MMP profile), show a similar distribution for both 
types of source measurements with most values between 0.1 and 
2. Therefore, the average UV radiation field energy density of the 
sources seems to be typically from a few tenth up to a factor 2 the 
intensity of the local Milky Way interstellar radiation field, as de- 
scribed by the MMP curve. 

The bottom right panel shows the distribution of Xcoi- As ex- 
plained in §5, Xcoi expresses the ratio of the optical to UV radi- 
ation field energy density (relative to the standard MMP curve). 
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Figure 5. Box-and-whisker plots showing the distribution of values for the compact source dust mass (top left), temperature (top middle), F24 (top left)), 
XUV (bottom right) and Xcol (bottom middle). In each panel the source measurements with high background (left) and low background (right) subtractions 
are shown. The solid horizontal line denotes the median value with the box delineating the 25th and 75th quartile range. The lines extend to the maxima and 
minima with open circles marking outliers (see text). 



Here the two measurements show some differences with the HB 
source measurement showing higher values and a greater range. 
This could mean that the average optical/UV ratio needed to fit 
the source emission is higher when one subtracts the higher GET- 
SOURCES background. Alternatively, this could mean that the PAH 
abundance required for the HB fit is lower than for the LB fit (be- 
cause of the Xcoi-PAH abundance degeneracy described in §5). In 
other words, a lower fraction of 8 /^m emission relative to the FIR 
emission is removed by the background estimated in the interarm 
apertures. 

As explained in §5, the cold dust temperature is not fit in 
our SED fitting procedure but is determined by the values of both 
Xuv and xjcoi. Both source measurements show similar tempera- 
ture ranges (see top middle panel of Fig. [5}, with the HB measure- 
ment only slightly higher (about 1 K higher). 

The source dust masses are shown in the upper left panel. In 
both cases, the compact sources have a quite narrow range in dust 
masses peaking close to 10° Mq and with most values between 
10 4 and 1O 6 M0. The small range of masses is likely due to our 
resolving power, which prevents us from detecting fainter sources 
which tend to also have smaller dust masses. The LB measurements 
for the sources have higher masses, because the background deter- 
mined in the interarm aperture subtracts less flux than the GET- 
SOURCES background. 

The distribution of F24 , the relative amount of 24 /im emission 
associated with PDR, is shown in the top right panel. Both source 
measurements show relatively similar values, mainly in the range 
0.25-0.6, which suggests that a substantial part of the 24 /im emis- 
sion is contributed by stochastically heated small grains in the dif- 
fuse dust component. However, we note that PDR SED dust emis- 



sion spectra can vary substantially (Groves et al. 2008) and might 
not be accurately reproduced by our PDR template for each indi- 
vidual case. 



7.7.2 Source Luminosities 

The SED fitting procedure also provides us with a measure of the 
total infrared luminosity of the sources as well as the luminosities 
due to the PDR and the diffuse dust emission component. Fig. [6] 
shows the distribution of the inferred luminosity values for each 
component. The total luminosity values are found in a rather small 
interval around 10 41 ergs/s. This is likely due to the detection tech- 
nique. We find that GETSOURCES does not detect sources below a 
minimum dust infrared luminosity of fa 0.5 x 10 41 ergs/s. Since 
there are not many sources with dust luminosities higher than a few 
times 10 41 ergs/s, the inferred dust luminosity range covered is 
rather small and this has important consequences on the inferred 
range of SED fitted parameters and SFR values (see discussion in 
§8). Fig. [6] also shows that the diffuse dust emission component 
dominates the dust luminosity. We point out that this does not nec- 
essarily mean that the dust heating from star forming regions is not 
responsible for most of the dust emission, since the diffuse dust 
emission component can be powered by both radiation from older 
stellar populations and the fraction of UV photons escaping from 
PDRs. The origin of the dust heating will be further discussed in 
§8.3. 
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Figure 6. Box-and-whisker plots showing the distribution of luminosities 
for the sources in the high background (left) and low background (right) 
subtractions. We show the distributions of the total luminosity, PDR lumi- 
nosity component and total diffuse component. See Fig.[3]for details. 



7.1.3 Star formation rates and efficiencies 

Fig. [7] shows the distribution of the SFR and SFE inferred val- 
ues. Due to the small areas considered here (i.e. average radius of 
sources is 150 pc), the SFRs are quite low. As outlined in §6 if 
the SFRs are too low (i.e. below 0.001 Mq yr~ ) then the calibra- 
tion of Ha and 24 /im fluxes may not adequately trace the SFR. 
However, the majority of the sources have SFRs above this value 
and furthermore, we have restricted ourselves to regions associated 
with recent star formation. These regions should have a good cor- 
respondence between gas and star formation. We note that the Ha 
contribution represents half or slightly more than half of the SFR 
for the sources. The median ratio of the Ha luminosity, L(H a ), to 
the scaled 24 /im luminosity, 0.031Z/(24/im), is 1.4. 

The SFEs of the source measurements have median values 
with median deviation values of 3.57 ± 3.45 x 10~ 10 yr _1 and 
3.1 ± 3.1 x 10~ 10 yr _1 in the case of the HB and LB measure- 
ment, respectively. This is consistent with the findings of Leroy et 
al. (2008) who found that in terms of molecular gas, the SFE has a 
mean value of 5.25 ± 2.5 x 10~ 10 yr" 1 . 

The LB and HB measurements show considerable differences 
between the inferred SFR and source gas masses, with the LB mea- 
surements having higher values in both cases. However, the LB and 
HB measurements show roughly the same range of SFE values. 
This means that roughly the same relative amount of background is 
being subtracted for both the SFR tracers and dust emission. This 
suggests the amount of background subtraction does not affect the 
results of the SFE, provided the same technique is used for both the 
gas and SFR tracers. 
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Figure 7. Box-and-whisker plots showing the distribution of values for the 
SFR (left) and SFE (right). In each panel the source measurements with 
high background (left) and low background (right) subtractions are shown. 
See Fig.|5]for details. 
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Figure 8. The radial variation of the luminosity components including the 
total (black circles), PDR component (blue triangles) and diffuse compo- 
nent (red squares), for the compact sources with high background subtrac- 
tion (left) and low background subtraction (right). In the upper left corner 
we list the Spearman rank correlation coefficient for each component and 
radius. The relative fraction of the PDR component to the total (blue tri- 
angles) and diffuse component to the total (red squares) are shown on the 
bottom panel for both cases. 



In order to determine the deprojected radial position, we as- 
sume an inclination and position angle for M83 of 24° and 225° 
respectively (Tilanus & Allen 1993). 



7.2 Radial Variations 

We now turn to the radial variations of the SED fitting parame- 
ters, the SFR and SFE and of the dust luminosities. Previous works 
which have used azimuthal averages (e.g. Munoz-Mateos et al. 
2009, Pohlen et al. 2010, Engelbracht et al. 2010, Boquien et al. 
2011) have found that the parameters describing the dust, includ- 
ing temperature and surface mass density, and the SFR tend to de- 
crease with radius particularly for late-type spirals. These studies 
have averaged over both the compact regions and the more diffuse 
emission. By separating the compact regions, we can better disen- 
tangle how and if the radial position affects the dust properties of 
these sources. 



7.2.1 Source infrared luminosities 

Fig. H] shows how the luminosity of the sources (upper panels) and 
the relative fraction of PDR and diffuse dust emission luminosity 
(lower panels) varies with radius. In the upper left comer we list the 
Spearman rank correlation coefficient for the correlation between 
the luminosity components and radius. Values of the correlation 
coefficient approaching +1 or —1 reflect stronger correlations or 
anti-correlations respectively. Values close or less than 0.5 repre- 
sent mild or weak correlations. In both the HB and LB measure- 
ments, we find only a mild correlation between the source lumi- 
nosities and radius. The source luminosities show a slight decrease 
with radial position with scatter. We also note that the regions with 
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Figure 9. The radial variation of the mass, SFR and SFE for the two com- 
pact source measures. The upper right hand corner lists the Spearman rank 
correlation coefficient. The error bars in red denote the median uncertainties 
of the points in three bins between 0-2 kpc, 2-4 kpc and 4-6 kpc. 
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Figure 10. The radial variation of Xcob Xuv, -F24 and temperature for the 
two compact source measures. The upper right hand corner lists the Spear- 
man rank correlation coefficient. The error bars in red denote the median 
uncertainties of the individual measurements in three bins between 0-2 kpc, 
2-4 kpc and 4-6 kpc. 



the highest PDR luminosity fraction are found in the inner regions, 
at a radial distance of 2-3 kpc from the galaxy centre (roughly the 
end of the bar). However, the diffuse dust component dominates the 
emission for all the sources, except in a few cases. 



7.2.2 Radial Variations of the SED Fitting Parameters, SFR and 
SFE 

Fig. [9] and [10] show the SED fitting parameters and the SFR and 
SFE values plotted versus the source deprojected radial positions 
for both the flux measurements of the compact sources. The median 
uncertainties of the individual measurements are shown by the red 
error bars and the Spearman rank coefficient is listed in the upper 
right of each panel. 

For both the HB and LB measurements, we find little to no 
variations with radius for any of the SED parameters including dust 
mass, temperature, F24, Xuv, Xcol, There is a slight trend showing 
that the SFR decreases with radius, but this is far less than what 
is typically seen in studies on late-type spiral galaxies that employ 
azimuthal averages. It seems that the properties of the compact re- 
gions are quite uniform and do not vary much with location in the 
galaxy. This is perhaps not too surprising, since our source sample 
covers only a relatively small range of dust luminosities, as shown 
in §7.1. Therefore, the sources we consider likely represent similar 
types of objects. Although the low spatial resolution did not allow 
us to detect a larger range of source luminosities, it is striking that 
the detectable bright sources lie almost exclusively on the spiral 
arms or at the ends of the bar in M83, which means that the local 
environment may in fact be quite similar as well. This suggests that 
azimuthal averages of dust properties could well mask important 
differences in environment within an azimuthal bin. The radial de- 
creases that are typically seen in other galaxies and for M83 in F12 
using azimuthal averages might simply reflect the increasing con- 
tribution of the interarm regions at larger radii. We note, however, 
that the SFE is quite constant with radius. This finding is consistent 
with studies that use azimuthal averages and pixel-by-pixel anal- 
yses (e.g. Leroy et al. 2008, Blanc et al. 2009, Foyle et al. 2010, 
Bigieletal. 2011). 



7.3 Correlations between source properties 

Having seen how the dust and star formation parameters of the 
compact sources vary with radius, we now turn to how they are 
inter-related and what types of correlations may exist between 
some of these parameters. In particular, we examine how the source 
dust mass and temperature are related to the radiation field param- 
eters xuv and Xcol and with the star formation rate and efficiency. 
Fig. QT] shows the dust temperature, Xuv, SFR and SFE plotted 
versus the dust mass, while Fig. 1121 shows xuv, X^oi, SFR and 
SFE values plotted against the dust temperature^. As before, the 
red bars show the median uncertainties in equally spaced bins and 
the Spearman correlation coefficients are shown in the upper-left of 
each panel. In general, we find that both LB and HB measurements 
of the compact sources present quite similar results. 

We observe a strong correlation of xuv with the dust temper- 
ature, while the optical-UV intensity ratio Xcol is not as strongly 
correlated to the dust temperature measured for the sources . These 
findings reflect the fact the UV radiation field intensity is the fun- 
damental parameter that determines the average dust temperature, 
while the optical-UV intensity ratio varies randomly among sources 
having same average dust temperatures. We point out that this does 
not mean that optical photons are not an important source of dust 
heating. Actually, the optical part of the radiation field accounts for 



6 In Appendix C we present how these parameters vary as a function of 
colour. 
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The Spearman rank coefficient is shown in the upper right of each panel. 
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more than 50% of the dust heating for the diffuse dust SED com- 
ponent when Xcoi > 1. 

The upper left panel of Fig.llllreveals a mild anti-correlation 
between the dust temperature and mass for the sources. The typ- 
ical error bars are much smaller than the inferred range of dust 
mass and dust temperature. Thus, this anti-correlation cannot be 
explained by the uncertain determination of dust temperature and 
mass, which are connected in the SED fitting since approximately 
Ldust oc MdustT^Jt . Therefore, we are confident that the sources 
with higher dust masses tend to have lower dust temperatures in the 
sample we considered. Given the strong correlation between \uv 
and temperature, we also find an anti-correlation between xuv and 
dust mass in the second row from the top of Fig.[TT] We will discuss 
the significance of the dust mass - dust temperature anti-correlation 
in §8.3. 

Fig. QT] (third row, left) shows that for both measures of the 
compact sources the dust mass is only weakly correlated with the 
SFR. If the dust traces gas, this plot should be similar to the S-K re- 
lation (Kennicutt 1998), except that typically surface densities are 
plotted rather than total values. The S-K relation for nearby galax- 
ies typically shows a tight correlation between gas surface densities 
and the SFR surface densities, at least when the scales considered 
are large enough (e.g. Schruba et al. 2010, Bigiel et al. 2008, see 
other references in §2.1). The relatively constant SFR with mass 
suggests that more massive sources are less efficient at forming 
stars and will be further discuss in §8.2. We also find that for both 
measures of the compact sources, the SFR is only mildly corre- 
lated with the dust temperature (see third row of Fig. 112) . This is 
not unexpected since, as explained in Sect. §2.2, the SFR is con- 



nected to the dust temperature through the dust mass and the total 
luminosity- weighted optical-depth. 

In contrast to the lack of correlations for the SFR, we find a 
mild and strong correlation respectively of dust mass and dust tem- 
perature with the SFE. While the SFE was roughly constant with 
radius, we find that the SFE is anti-correlated with dust mass and 
correlated with dust temperature. These findings have not been seen 
before, as most studies have found that the SFE does not vary much 
with other properties in terms of pixel-by-pixel and azimuthal aver- 
ages (i.e. Leroy et al. 2008). This will be further discussed in §8.2 
and §8.3. We note that that the anti-correlation between the SFE 
and mass is primarily due to the fact that the SFR does not vary 
much with the inferred gas mass. Thus, the correlated axes produce 
the anti-correlation. To check this, we performed Monte-Carlo sim- 
ulation assuming the SFR was roughly constant with a Gaussian 
distribution and a mean and standard deviation equivalent to that 
found in our measurements. We assumed a similar distribution for 
the simulated source masses. The simulation naturally produced a 
linear plot of SFE versus mass with a slope of -1.05 ± 0.15. The 
compact sources in our analysis produced a relation with a slope of 
-0.71 ± 0.14. The slopes while not identical within the uncertain- 
ties, are similar, suggesting that the correlated axes are the primary 
reason for the relation. 

We also examined the relation between dust luminosity and 
the SFR, as shown in Fig[T3] We see that not only is the total dust 
luminosity but also the PDR and diffuse dust emission luminosities 
are well-correlated with the SFR. As explained in section §2.2, this 
is what is expected if a young stellar population is responsible for 
the dust heating and the total luminosity-weighted opacity of the 
sources, r, does not vary arbitrarily among the sources. 
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Figure 13. Dust luminosity versus the SFR for both measures of compact 
sources including the total luminosity (black circles), PDR component (blue 
triangles) and diffuse component (red squares). The correlation coefficients 
for each component are shown in the upper right. The solid line denotes the 
relation of SFR=2.8xlO~ 44 x L dust (Calzetti 2012), which assumes that 
all the young stellar luminosity is absorbed and re-emitted by dust. 



8 DISCUSSION 

In this section we discuss the properties of the detected compact 
sources of M83 in terms of gas masses and locations (§8.1), star 
formation (§8.2) and dust heating (§8.3). 



8.1 Source Gas Masses and Locations 

The inferred dust masses of the compact sources lie within the 
range of 10 4 -10 6 Mq. These dust masses can be used to estimate 
the source gas masses as discussed in §6. We use a gas-to-dust ratio 
of 100, which means the gas masses of the sources have values in 
the range 10 6 -10 8 M©. These values correspond to the high end 
of the molecular cloud mass distribution function derived both ob- 
servationally (see e.g. Solomon et al. 1987, Rosolowski et al. 2005, 
Gratier et al. 2012) and by numerical simulations (e.g. Nimori et al. 
2012). Furthermore, previous CO studies of nearby galaxies have 
highlighted the presence of very massive clouds with gas masses 
of order of 10 7 Mq or higher, which are usually referred to as Gi- 
ant Molecular Associations (GMA) (e.g. Vogel et al. 1988, Koda 
et al. 2009, Muraoka et al. 2009). The high masses of our com- 
pact sources suggests that they are GMAs rather than less massive 
GMCs. 

However, because dust continuum emission does not provide 
kinematical information, it is not possible to check whether the de- 
tected sources are gravitationally bound clouds. The smallest cloud 
radius we can measure is 130 pc due the 70 fim beam size (6"). The 
majority of the detected sources are found on the spiral arms. Stud- 
ies of the molecular gas in M83 have shown that most of the bright 
compact sources on the spiral arms are bound clouds in GMAs. 
Muraoka et al. (2009) examined CO in M83 and measured the virial 
parameter, a, defined as the ratio of the virial mass to the CO lumi- 
nosity mass. They found that a is almost equal to unity in on-arm 
clouds, suggesting that they are in a gravitationally bound state (see 
also Rand et al. 1999, Lundgren et al. 2004). The resolution of their 
study was 7". 5, which is comparable to our resolution. Thus, it is 
quite likely that our sources are bound GMAs. 

The spiral arms are the natural location where ISM gas is 
brought to higher densities, thus leading to the formation of mas- 
sive clouds. Therefore, it is not surprising, given our resolution and 
sensitivity, that we only detect compact sources on the spiral struc- 
ture. We note that the bar-spiral arm transition region harboured the 
most sources and that just few were detected along the length of the 
bar. Due to shear motions in bars, GMCs can be easily pulled apart 



(Downes et al. 1996). In contrast to spiral arms, bars exhibit lower 
star formation rates and in some cases lower SFE (Momose et al. 
2010). Meanwhile, studies of atomic and molecular gas have shown 
that the most massive complexes are found in the bar-spiral arm 
transition region. This can be explained by orbit crowding when gas 
on highly elliptical orbits in the bar converges with the gas orbit- 
ing in the spiral structure (e.g. Kenney & Lord 1991). The location 
of the compact sources detected in the FIR is consistent with these 
findings in the gas. 

8.2 Star Formation Rates and Efficiencies 

One of the goals of this work has been to infer cloud gas masses 
using the dust emission of compact sources. Combined with a mea- 
sure of the SFR, we can then examine the SFE in these regions. 

There have been many recent studies that have examined star 
formation on spatially-resolved scales. Typically these studies have 
employed one of three methods: pixel-by-pixel analyses, where the 
SFR tracers and gas tracers are compared in individual pixels; az- 
imuthal average analyses, where average values in radial bins are 
compared; or aperture photometry on compact sources bright in 
some SFR tracer such as Ha emission. 

Our work differs in several key ways and, before discussing 
the results, we summarize these differences here: 1) We use the 
dust emission to detect the compact sources. While most studies 
have focused on aperture photometry of active star forming regions 
(e.g. Calzetti et al. 2005, Kennicutt et al. 2007, etc.), few have se- 
lected compact sources as seen in dust emission maps, which may 
be more akin to locating peaks in the gas (e.g. Schruba et al. 2010). 
However, there is a good correspondence between the detected FIR 
sources and star formation regions, which is demonstrated by the 
tight correlation between total infrared luminosities of the sources 
and their SFRs (see Fig. 13); 2) we perform a diffuse background 
subtraction of the dust emission using both a local background de- 
fined by GETSOURCES and one determined from nearby apertures 
in the interarm regions. Previous studies have performed a back- 
ground subtraction on the star formation tracers but not on the gas 
mass (though see Rahman et al. 201 1); 3) it is important to keep in 
mind, when comparing this work with others, that we analysed the 
total SFR and dust masses for individual sources while most stud- 
ies consider the surface densities of the gas mass and SFR when 
examining trends between the two. The reason why we did not con- 
sider surface densities is that it is not clear what the source area 
should be, since the apparent size of the sources varies on each 
map depending on the map resolution. In addition, since most of 
the sources are not well-resolved, it is not possible to determine the 
physical source area which is typically smaller than the PSF beam 
area. 

As discussed in §2, previous works have found that, on rela- 
tively large spatial scales (> 500 pc) clear correlations are found 
between the gas, SFR and dust as well as a tendency for these 
parameters to decline with radius. However, as we saw in §7, the 
compact regions detected in this study show a relatively constant 
SFR and dust mass with radius, albeit with large scatter and there 
is little to no correlation between the SFR and dust mass. Due to 
the small spatial scales of these sources (< 300 pc), one does not 
expect to recover the relations found by averaging over much larger 
areas. Recently, it has become clear that on small spatial scales (less 
than 400 pc ) the S-K relation between the gas surface density and 
SFR surface density breaks down (Calzetti et al. 2012, Feldmann 
& Gnedin 2011 and Schruba et al. 2010). This has been found to 
be even more prominent in the case where gas peaks are selected 
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(Schruba et al. 2010). The scatter is due to the fact that individ- 
ual GMCs and HII regions will show varying gas mass to star for- 
mation rates depending on their evolutionary state. Averaging over 
large areas means that many objects in different states are averaged 
and the relation between SFR and gas mass is recovered. 

The scales on which the SFR and gas mass surface density are 
measured play an important role in defining the relation (Liu et al. 
2011). Typical SFR tracers suitable for relatively large galaxy re- 
gions are not straightforwardly applicable to small regions of sub- 
kpc scales. In fact, in this case, the intrinsic assumption that the 
stochastic characteristics of star formation are averaged out from 
the integration on large areas can easily break down (i.e. Calzetti et 
al. 2012). This happens because, one ideally needs a complete sam- 
pling of the stellar IMF. Small regions may not encompass a large 
enough stellar population to do so. Furthermore, beyond the issue 
of sampling the stellar IMF, there is the problem of time-averaging. 
The SFR is calibrated based on the assumption of a constant SFR 
over 100 Myr (Calzetti et al. 2007). While this is true for entire 
galaxies, on small scales this assumption may no longer hold and 
can create a large scatter in the relation. 

Our sources have sizes of roughly 300 pc and thus, we find 
large scatter between the inferred gas mass and the SFR. We also 
note that our sources populate only a narrow range in FIR lumi- 
nosity and dust mass (see Fig|5]and[6]l. Thus, we see a relatively 
constant SFR for sources with such similar characteristics. If we 
could populate the plot with lower mass regions, an S-K relation 
with considerable scatter might be recovered. Our results reflect 
stochasticity and the fact that individual star forming regions can 
exhibit a range of SFRs for a given dust or gas mass. However, we 
should note a few other caveats that might introduce scatter into the 
relation. 

First, it is possible that dust emission does a poor job of tracing 
molecular gas especially in light of our use of a constant GDR. 
There may be large variations in the GDR (i.e. due to metallicity 
variations), which means that our gas mass estimates might not be 
reliable. However, the compact regions all lie on the spiral arms 
and beyond the nuclear region of the galaxy. In F12, the largest 
variations in the GDR were seen in the central regions and, beyond, 
the GDR was quite constant. It is also possible that some of the 
peaks found in the dust emission are not associated with active star 
forming regions. However, we found a correlation between the dust 
luminosity and the SFR and studies like Verley et al. (2010) have 
found a tight correlation between compact regions in M33 and star 
formation rate tracers. We should also note that all of the sources 
were detected in the 70/im map, which is also known as a good 
tracer of the SFR (Li et al. 2010). 

A second possibility is that our background emission subtrac- 
tion has introduced scatter. Rahman et al. (2011) showed that as 
one increases the amount of diffuse emission subtracted from the 
gas, the scatter in the S-K relation increases. However, particularly 
given that we are using dust emission, which may trace gas of dif- 
ferent phases, it is important to account for a diffuse component. 

Despite the large scatter, we find that the SFR varies little with 
the inferred gas mass and is roughly constant. This implies that the 
more massive clouds are less efficient at forming stars. A physical 
reason for this is that the more massive clumps might be more ex- 
tended and thus have lower gas surface densities. Studies of GMCs 
have revealed that only the densest gas is directly associated with 
star formation. In this way, the more massive and thus extended 
regions may have roughly comparable SFRs (or lower SFEs). We 
note that we also found that there is a strong anti-correlation be- 
tween the SFE and the inferred gas mass. However, as discussed in 



§7.3, this is is largely a product of the correlated axis since mass 
appears both in the y and x-axes (see Murray (201 1) who examined 
galactic GMCs and found an anti-correlation which was attributed 
to a similar effect). We also note that we find no radial variation of 
the SFE of the compact regions. Previous works, which have exam- 
ined kpc-sized regions, have found that the SFE is relatively con- 
stant regardless of the variable considered (i.e. Leroy et al. 2008). 
This is consistent with our findings with radius. 



8.3 Dust Heating 

We showed in §2.2 that the simultaneous presence of a correlation 
between the SFR and source dust luminosity as well as a corre- 
lation between the SFE and the dust temperature suggests that the 
radiation impinging on the dust mass is mainly produced by a local, 
young stellar population. Contrary to common lore, a strong corre- 
lation between SFR and dust temperature should not be necessarily 
expected in this case. Indeed we saw in §7 (see Fig.[T3]and[T2]) that 
while the dust temperature was only mildly correlated with SFR, 
it was strongly correlated with the SFE. We also saw that the dust 
luminosity was correlated with the SFR (see Fig.ll3t. This supports 
a scenario where a young, local stellar population is powering the 
dust emission of the FIR bright sources of M83. 

In Fig. [T3] the solid line illustrates the maximum level of dust 
luminosity that can be powered by the young stellar population by 
assuming that all the luminosity from this population is absorbed 
and re-emitted by dust (optically thick case) using the following 
relation by Calzetti (2012): 



SFR = 2.8 x 10" 44 L d 



(6) 



This relation has been derived using the same assumptions on the 
star formation history, metallicity and IMF as the Ha based cali- 
bration we used to infer SFR of the sources. 

This line should delineate the maximum value, yet we find 
many sources lie close to or above this limit. Specifically we find 
that the percentage of sources that lie above or only within 2a be- 
low this limit is 60% or 20% in the HB or LB case respectively. 
While it is surprising to find sources that lie above or close to this 
maximum limit, these findings can be explained in one of two ways. 

Firstly, it is possible that some regions are indeed very close 
to being completely optically thick and therefore the total dust lu- 
minosity will be approximately equal to the intrinsic young stellar 
population luminosity in those cases. However, this is unlikely be- 
cause these regions are also Ho emitters, which means, by defini- 
tion, that at least some fraction of the young stellar population lu- 
minosity is able to escape unabsorbed (at least through re-emission 
in Ha). Furthermore, we checked that the differences between the 
observed dust luminosities and the dust luminosities for an opti- 
cally thick case, predicted by the relation in Eq. [6] do not depend 
on the ratio F(Hq)/SFR[Hq — 24/xm] . This suggests that the scatter 
of the observed points is not driven by differences in attenuation. 

Secondly, it is possible that there is a local older stellar popu- 
lation or stellar populations outside the projected source area con- 
tributing substantially to the heating. In this case, the approxi- 
mate equality between dust and young stellar population luminosity 
would be reached through this extra-heating from other radiation 
sources. Thus, the infrared total luminosity could still be used as a 
star formation rate indicator in those cases despite the fact that the 
dust heating is not powered exclusively by a local, young stellar 
population (see Bendo et al. 2012 also). This might explain why 
the recent work of Verley et al. (2010) and Boquien et al. (2010) 
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showed that the single Herschel 100, 160 and 250 /im band lumi- 
nosities of individual sources in M33 are good tracers of SFR, since 
most of the dust luminosity is emitted at FIR wavelengths. 

We should note that the SFE vs Tdust correlation is at least 
partially driven by the Mdust vs Td us t anti-correlation (more evi- 
dent for the LB measurement), since the SFR alone does not corre- 
late with A/dust and only mildly with Td us t . There are multiple po- 
tential explanations for an anti-correlation between mass and tem- 
perature. One possible explanation is that it is due to the different 
evolutionary stages of the sources, with the more massive sources 
having a relatively lower number of young, recently formed stars. 
In this way, the radiation field would be less intense for massive 
sources and the dust would be cooler. Alternatively, it could be that 
the more massive sources are more efficient at shielding radiation 
originating from outside the local region. As outlined in the pre- 
vious paragraph, there is likely some heating due to other stellar 
populations other than just the local, young population, and thus 
such shielding could keep more massive sources cooler. 

The anti-correlation between dust temperature and mass could 
be in part due to the lower limit in dust luminosity for the sources. 
FigQl] shows the inferred dust masses versus dust temperatures 
with the plotted points colour coded depending on the source total 
dust luminosity. Since most of the dust luminosity is emitted at FIR 
wavelengths where Ldust oc -MdustTdust , the points correspond to 
higher total dust luminosities while moving from the bottom-left to 
the upper-right region of the plot. From this figure one can see that 
there is a lower limit in luminosity (~ 5 — 6 x 10 40 erg/s), deter- 
mined by the minimum fluxes for source detection. This explains 
why there are no sources populating the lower-left part of the dia- 
gram. The steepness of the observed anti-correlation is consistent 
with the 4 + P exponent (with /?=2), as shown by the line in Fig. 1141 
However, higher mass sources do not populate uniformly the en- 
tire range of temperatures in the upper-right part of the diagram, 
where source detection would easily be possible. Therefore, while 
the low-luminosity limit may contribute to the anti-correlation be- 
tween dust temperature and mass, it cannot fully account for it. 
Thus, there should be at least some physical restrictions at play. 

We note that the lower limit to the luminosity range explored 
in our analysis does introduce a bias that should be taken into 
account. Essentially, our sources are the brightest FIR compact 
sources in the disk of M83 and they are associated with the high- 
end of the molecular cloud mass distribution. As discussed in §8.1, 
the sources are associated with gas masses in the range ~10 6 - 
10 8 Mq, while their total dust luminosities are in the range ~0.5- 
5x 10 41 erg/s. The results of this work should be considered valid 
only within these parameter ranges. 

We found a mild correlation between SFR and Td ust . How- 
ever, as shown in Appendix C, a better correlation is found when 
plotting the inferred 70/160 fim flux ratio versus the SFR. The 
difference between the two cases is due to the fact that Td ust in 
our SED fitting procedure refers to the cold dust temperature of the 
diffuse dust component. Meanwhile the observed 70 and 160/iin 
fluxes have contributions from both the diffuse component and the 
PDR component. Nonetheless, in Appendix C we also show that 
the correlation between SFE and 70/160 /im flux ratio is clearly 
tighter and this is consistent with the arguments presented in §2.2 
and very similar to the strong correlation found by plotting SFE 
versus T duEt . 

Pixel-by-pixel analyses have tended to find varying results 
with no strong consensus for a correlation between SFR and modi- 
fied black body dust temperatures or color temperatures. We should 
note that the pixel-by-pixel analyses have considered the star for- 




Figure 14. Source dust masses versus dust temperature for the HB (left) 
and LB (right) measurements. The points are colour-coded based on their 
total dust luminosity. The solid line represents the relation L dust °c Mdust 
T dust 6 for L dust =6x 10 40 ergs/s. 



mation rate surface density (hereafter SFRSD) as opposed to the 
SFR used here. For M83, F12 found a correlation between the 
SFRSD and Td ust . A similar result for M83 has been found by 
Bendo et al. (2012) when comparing the 70/160 /im flux ratio with 
the observed Hq brightness, uncorrected for internal dust attenua- 
tion. However, Bendo et al. (2012) found a weaker correlation for 
M81 and NGC 2403. H Boquien et al. (2011) found a mild corre- 
lation between SFRSD and the 70/160/xm flux ratio for M33. In 
addition, Smith et al. (2012) for the Andromeda galaxy and Skibba 
et al. (2012) for the Large and Small Magellanic Clouds found no 
or only a weak correlation between SFRSD and Td ust . Therefore, it 
seems that in general the correlation between SFR vs tracers of dust 
temperature is not observed to be very strong when one considers 
either total SFRs for individual FIR sources or SFR averages on ar- 
eas equivalent to pixel sizes. As explained in §2.2, this is expected 
even for the case where young stellar populations are responsible 
for the dust heating because the dust temperature also depends on 
the dust mass and the total luminosity-weighted optical depth (see 
Eq.[4j), which are not uniform throughout a galaxy disk. 

As we have seen, while isolating FIR bright sources has high- 
lighted the presence of interesting correlations between the inferred 
source parameters, the narrow range of luminosities of the sources 
has introduced a potentially important bias. In order to elucidate 
the influence of this bias on the inferred relationships, we need to 
expand this work to include sources within a wider range of lumi- 
nosities, particularly fainter sources. To do so, we plan a follow-up 
study which will include other spiral galaxies that are even closer 
(i.e. NGC 2403 and M33). The closer proximity will allow for an 
even better spatial resolution, which is necessary to detect sources 
in a wider range of luminosities and determine if the results for the 
GMAs in M83 are also valid for GMCs in other galaxies. 



9 CONCLUSIONS 

The main aim of this paper has been the investigation of the star 
formation and dust heating properties of the compact FIR bright 
sources as observed on the Herschel maps of the nearby spiral 
galaxy M83. By combining the source detection and photometry 
algorithm GETSOURCES, the dust emission SED fitting method 
of Natale et al. (2010) and the [Ha - 24 /im] star formation rate 



7 The weaker correlation is attributed to the dust heating by bulge stars in 
M81 and artefacts in the images of NGC 2403; see the discussion in that 
paper. 
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calibration by Calzetti et al. (2007), we have developed a new 
procedure to determine gas masses, radiation field intensities, cold 
dust temperatures, dust luminosities, star formation rates and star 
formation efficiencies associated with those sources. 

The main results of our analysis are the following: 

• We have found that the well-detected compact FIR sources 
are mostly associated with giant molecular associations, with gas 
masses in the range 10 6 -10 8 M© and dust total infrared luminosi- 
ties in the range 0.5 x 10 41 -10 42 ergs/s. The majority of the sources 
are located on the spiral arms of M83 with only few sources found 
in the interarm and within the bar in the central region of M83. 

• None of the inferred physical quantities for the sources shows 
a strong variation with radius, including SFR, gas masses and dust 
temperature. Previous studies have usually found radial variation 
for these quantities, although only after averaging on larger areas 
including interarm regions and without subtracting any local back- 
ground. 

• The SFR does not seem to correlate strongly with the gas mass 
of the sources. The lack of correlation is most likely due to the small 
spatial scales considered (~200-300 pc) and/or the relatively small 
range of inferred gas masses. 

• The star formation efficiency SFE, defined as SFR/M gas , 
shows an anti-correlation with source gas mass. This finding sug- 
gests that the more massive GMAs are less efficient in forming stars 
in the last few Myr. However, we note that this anti-correlation is a 
consequence of the roughly constant SFR with inferred gas mass. 

• We found that the SFR correlates well with total dust luminos- 
ity, which is consistent with a scenario where dust is predominantly 
heated by the local young stellar population. However, between 20- 
60% of the sources show dust luminosities which are greater than 
or only 2a below those predicted when the heating is due only to a 
local young stellar population, embedded in a optically thick dust 
distribution. It is unlikely that the sources are completely optically 
thick, as they are also bright in Ha. Thus, it seems that there must 
be some extra-heating by either a local older stellar population or 
by an external radiation field. 

• We found a correlation between the SFE and Td us t which is 
tighter than the mild correlation we found between SFR and Td ust . 
This is expected if the dust is heated primarily by recent star for- 
mation. 

• We found a mild anti-correlation between dust and mass tem- 
perature. While our sources have a low-luminosity limit, which 
may contribute to this anti-correlation, we find that this can not 
fully account for it. Thus, we speculate that the more massive 
sources are more efficient at shielding from an impinging radia- 
tion field or that more massive sources are in an earlier stage of star 
formation. 

We plan to use the same procedure presented in this pilot 
work on a set of nearby galaxies observed by Herschel. This will 
help to further clarify the origin of the observed correlations and 
their implications for star formation and dust heating associated 
with FIR bright sources on galactic scales. 
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APPENDIX A: RESOLUTION TESTS 

We test for any systematic errors in the process of detecting and 
measuring fluxes of the compact sources using GETSOURCES. To 
do so, we take the 70 /im map with 6" resolution and degrade its 
resolution to match that of 160 /im map (12"), the 250 /im map 
(18.2") and the 350 /im (24.5") using the convolution kernels of 
Aniano et al. (2011). We then rerun the extraction process using 
these series of maps. We can then compare how the flux of the 
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Figure Al. Total flux of the sources measured by GETSOURCES using the 
degraded 70 /jm maps versus the total flux of those sources measured in 
the 70 fim map in its native resolution. The median uncertainties for each 
resolution are shown on the right. At poorer resolutions, the source flux is 
increasingly overestimated. 



sources is affected by the degradation of the resolution. In princi- 
ple, at each resolution, the flux should be the same for each source. 
However, particularly in crowded regions, the source flux is more 
complicated to determine as the resolution is degraded. Fig. IA1I 
compares how the source flux at the degraded resolution compares 
to the original. In general, we find that the flux is somewhat overes- 
timated as the resolution is degraded. Half of the sources fall within 
the uncertainties on the measurements and the majority lie within 
a factor of 2 of the flux determined at the 6" resolution. We also 
tested how such a systematic uncertainty could affect the results of 
the inferred dust temperature and mass. By applying the average 
correction factors to the source FIR flux measurements (see §4.2) 
and performing SED fitting as described in §5, we found an average 
relative difference in dust mass of ~ 50% and an average difference 
of 2K for the dust temperatures. Although non negligible, the pos- 
sible systematic effect highlighted by this test is not large enough 
to affect the main conclusions of this work. 



Figure Bl. Comparison of the dust temperature derived by the SED fitting 
to that obtained by a modified blackbody fit of the FIR wavelengths. The 
median uncertainties are shown on the bottom right. 
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Figure B2. Comparison of the dust mass derived by the SED fitting to that 
obtained by a modified blackbody fit of the FIR wavelengths. 



emission models of Draine & Li (2007) the total mass of M83 was 
greater by a factor of 1.3 than those found using a modified black- 
body fit. 



APPENDIX B: COMPARISON WITH MODIFIED 
BLACKBODY FITS 

We compare the results of our two component SED fitting method, 
by determining the dust temperature and mass using a modified 
blackbody fit to the FIR emission (70-350 /an). 
We use a modified blackbody of the form: 



Sx=NBi{T)\- 



(Bl) 



where S\ is the flux density, B^(T) is the Planck function, N is a 
constant related to the column density of the material and /? is the 
dust emissivity index. We hold f3 at a constant value of 2, based on 
the average value found by F12. We fit the points using MPFIT, a 
least-squares curve fitting routine for IDL (Markwardt 2009). We 
refer to the temperature and dust masses determined by this tech- 
nique as 1 dust, mb and M du st,MB respectively. 

Fig. IB II and IB2I compares the temperatures and masses de- 
rived from the SED fitting with those from a modified blackbody 
fit. There is close agreement with the temperatures, although the 
modified blackbody dust temperatures tend to be slightly higher. 
The masses found by the SED fitting are larger than those from 
the modified blackbody fits with a median ratio of a factor of 1.4. 
Dust emission models typically have higher masses than modified 
blackbody fits, because the latter do not encompass the full range 
of temperatures of the dust as they use only the FIR maps (Dale 
et al. 2012). In F12, they found that fitting the SED using the dust 



APPENDIX C: COLOR CORRELATIONS 

Many previous works that have examined dust heating mechanisms 
and correlations with recent star formation rate tracers have used 
dust emission colors as opposed to SED fitting parameters (i.e. 
Bendo et al. 2012, Boquien et al. 201 1). For the purposes of com- 
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Figure CI. Ratio of 70-to-160 fim emission (top), 160-to-250 /im emission 
(second from top) and 250-to-350 (im emission (bottom) for the compact 
source measurements with high background subtraction (left) and low back- 
ground subtraction (right) versus deprojected radial position. 
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Figure C2. Ratio of 70-to-160 (im emission (top), 160-to-250 /im emission 
(second from top) and 250-to-350 fim emission (bottom) for the compact 
sources with high background subtraction (left) and low background sub- 
traction (right) versus the star formation rate measured in each case. 



160 ^m intensity ratio. This correlation is stronger than that seen 
between the SFR and the dust temperature Td us t inferred by SED 
fitting. The difference in the results when using the 70-to- 1 60 /im 
intensity ratio and Tdust is due to the fact that the latter is the dust 
temperature associated with the cold diffuse dust component in the 
SED fitting. At both 70 and 160/xm, a part of the emission is con- 
tributed by the PDR component as well, thus the FIR colors fitted 
by the diffuse dust component are usually slightly different than the 
observed ones. 

However, we note that there is a strong correlation between 
the SFE and the 70-to-160 jim intensity ratio (see Fig, IC3t . which 
is very similar with what we saw in §7.3 for the correlation ob- 
served between SFE and Td u8 t- On the other hand, we did not find 
any strong correlation between 160/250 /jm or 250/350 fim flux ra- 
tios with SFR or SFE. This might indicate that the dust emission 
at these wavelengths is more affected by heating from older stellar 
populations. However, one should also note that the flux uncertain- 
ties are quite large at longer wavelengths and, as a consequence, 
there could also be a stronger scatter hiding the presence of intrin- 
sic correlations. 
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Figure C3. Ratio of 70-to-160 fim emission (top), 160-to-250 fim emis- 
sion (second from top) and 250-to-350 fim emission (bottom) for the com- 
pact source measurements from GETSOURCES (left) and alternative com- 
pact source measurements (right) versus the SFE measured in each case. 



APPENDIX D: SOURCE FLUXES 



parison, we also examine how the dust emission colors vary with 
radius and with the star formation parameters SFR and SFE. 

We consider the color ratios of 70-to-160 /im, 160-to-250 /im 
and 250-to-350 /jm emission. In following what we have done in 
the previous sections, we consider both measures of the compact 
sources. 

Fig. |C1| shows how the colors vary with radius. As we saw 
with the SED fitting parameters, for both measures of the com- 
pact sources, we find no or just a weak trend in the colors with ra- 
dius. Both Bendo et al. 2012 and Boquien et al. 201 1 found in their 
pixel-by-pixel analysis that the FIR colours present decreasing ra- 
dial trends. Our measurements show that, when compact regions 
are separated, they do not seem to behave in the same way. 

Fig. |C2| shows how the colours correlate with the SFR. Bendo 
et al. 2012 and Boquien et al. 2011 highlighted the presence of 
correlations between SFR and FIR colors, especially the 70-to-160 
/jm intensity ratio. This finding was used to imply that at shorter 
FIR wavelengths the dust heating is powered mostly by recent star 
formation. We find a mild correlation between the SFR and 70-to- 
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Table D 1 : HB measurement source fluxes and positions 





RA 


5 


r8/im 


F24/jm 


F7UMm 


Fl60M"J 


F250m»™ 


F350 M m 


F(H a ) 




[deg] 


[deg] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[MJy] 


1 


204.27867 


-29.826725 


79 ±6 


190 ± 10 


2710 ± 150 


3690 ± 410 


1640 ±410 


720 ± 220 


51 ±2 


2 


204.26433 


-29.900421 


35 ±5 


139 ±9 


1495 ± 92 


2490 ± 290 


1320 ±310 


650 ±210 


36 ±2 


3 


204.22128 


-29.858999 


89 ±8 


180 ± 10 


2410 ± 160 


3790 ± 680 


1830 ±480 


650 ± 220 


82 ±3 


4 


204.28271 


-29.854724 


68 ±7 


107 ±7 


2090 ± 130 


3060 ± 500 


1540 ± 450 


410 ± 200 


34 ±3 


5 


204.23187 


-29.830846 


88 ±7 


73 ±6 


1860 ± 120 


3120 ±370 


1440 ± 410 


580 ±210 


84 ±3 


6 


204.22652 


-29.884745 


70 ±6 


100 ± 10 


2010 ± 140 


2530 ± 600 


910 ± 360 


330 ± 180 


98 ±3 


7 


204.22894 


-29.885796 


48 ±6 


80 ±9 


1840 ± 140 


1550 ± 560 


390 ± 310 


- 


80 ±3 


8 


204.22172 


-29.882527 


36 ±5 


81 ±8 


1580 ± 120 


1900 ± 500 


1620 ± 430 


710 ± 220 


40 ±3 


9 


204.21930 


-29.863860 


120 ± 10 


140 ± 10 


2840 ± 200 


5690 ± 960 


1830 ±490 


1370 ± 290 


30 ±4 


10 


204.22046 


-29.880233 


102 ±8 


140 ± 10 


1630 ± 120 


3160 ± 610 


1890 ±460 


590 ± 210 


79 ±3 


11 


204.21971 


-29.852660 


53 ±6 


62 ±9 


1180 ± 110 


2110 ±340 


1230 ± 420 


460 ± 200 


65 ±4 


12 


204.26957 


-29.849437 


87 ±8 


106 ±9 


2260 ± 160 


3820 ± 730 


1900 ± 470 


800 ± 230 


127 ±3 


13 


204.28487 


-29.869891 


84 ±8 


71 ±6 


1427 ± 98 


2430 ± 380 


1280 ± 400 


500 ± 210 


63 ±3 


14 


204.22800 


-29.882334 


62 ±6 


58 ±8 


1620 ± 120 


3370 ±710 


2610 ± 530 


940 ± 240 


16 ±2 


15 


204.28634 


-29.858535 


31 ±5 


49 ±6 


531 ±60 


1970 ± 370 


980 ± 390 


370 ± 200 


8±3 


16 


204.28083 


-29.852275 


150 ± 10 


270 ± 20 


3850 ± 250 


5850 ± 840 


4590 ± 850 


1330 ± 300 


104 ±4 


17 


204.16814 


-29.855863 


29 ±2 


30 ±2 


554 ± 31 


940 ± 160 


490 ± 160 


216 ±74 


40.2 ± 0.9 


18 


204.22554 


-29.844806 


20 ±4 


60 ±5 


710 ±67 


1530 ± 290 


640 ± 370 


340 ± 190 


70 ±3 


19 


204.21463 


-29.883342 


25 ±5 


37 ±7 


666 ± 89 


770 ± 280 


1250 ± 350 


390 ± 190 


15 ±3 


20 


204.28397 


-29.881946 


37 ±6 


38 ±4 


820 ± 74 


1790 ± 280 


1010 ± 320 


430 ± 200 


42 ±3 


21 


204.20849 


-29.878634 


33 ±4 


42 ±5 


814 ± 92 


2420 ± 330 


960 ± 300 


530 ±210 


20 ±4 


22 


204.26997 


-29.823576 


50 ±6 


44 ±6 


986 ±91 


1830 ±320 


1030 ± 350 


300 ± 180 


86 ±3 


23 


204.29186 


-29.858028 


19 ±4 


12 ±4 


383 ±58 


580 ± 260 


730 ± 340 


220 ± 180 


43 ±4 


24 


204.25244 


-29.905010 


46 ±6 


23 ±3 


451 ±50 


1610 ±230 


1130 ±290 


570 ± 190 


38 ±2 


25 


204.21223 


-29.844608 


106 ±8 


99 ±8 


2040 ± 140 


3110 ± 410 


1390 ±410 


790 ± 230 


108 ±3 


26 


204.21877 


-29.863393 


16 ±4 


39 ±8 


478 ± 84 


760 ± 560 


1350 ± 420 


- 


83 ±3 


27 


204.28831 


-29.849784 


23 ±4 


25 ±4 


633 ± 65 


450 ± 260 


360 ±310 


- 


36 ±4 


28 


204.24629 


-29.907226 


34 ±3 


18 ±2 


453 ± 50 


1220 ± 210 


820 ± 250 


380 ± 170 


45 ±2 


29 


204.18835 


-29.879014 


16 ±2 


20 ±2 


350 ± 33 


810 ±200 


590 ± 210 


330 ± 110 


16 ±3 


30 


204.18091 


-29.872982 


12 ±2 


25 ±3 


317 ±39 


630 ±180 


540 ± 200 


119 ±76 


46 ±4 


31 


204.23124 


-29.845585 


69 ±7 


55 ±5 


730 ± 64 


1860 ±430 


590 ± 370 


310 ± 180 


18±2 


32 


204.29158 


-29.819528 


16 ±4 


17 ±2 


447 ± 61 


450 ± 250 


270 ± 140 


170 ± 110 


61 ±4 


33 


204.22354 


-29.813422 


13 ±2 


15 ±2 


300 ± 26 


540 ± 160 


360 ± 180 


340 ± 140 


44 ±2 


34 


204.29732 


-29.831076 


14 ±4 


16 ±2 


360 ± 47 


770 ± 280 


580 ± 230 


270 ± 150 


39 ±3 


35 


204.27795 


-29.823031 


12 ±3 


21 ±5 


241 ± 44 


1530 ± 300 


780 ± 290 


260 ± 170 


15 ±2 


36 


204.24819 


-29.810471 


59 ±5 


47 ±3 


908 ± 57 


1780 ± 240 


820 ± 260 


330 ± 150 


68 ±2 


37 


204.26738 


-29.899060 


13±4 


17 ±4 


381 ±46 


1590 ± 250 


790 ± 270 


410 ± 190 


28 ±2 


38 


204.23665 


-29.880081 


46 ±7 


52 ±6 


1040 ± 110 


2730 ±810 


1260 ± 380 


440 ± 190 


15 ±2 


39 


204.21163 


-29.866269 


23 ±5 


21 ±5 


423 ± 83 


460 ± 250 


430 ± 260 


150 ± 150 


20 ±4 


40 


204.31140 


-29.843049 


22 ±3 


14 ±2 


329 ± 33 


810 ±250 


400 ± 160 


250 ± 120 


- 


41 


204.25832 


-29.925232 


25 ±3 


16 ±2 


313 ±25 


790 ± 170 


410 ± 150 


230 ± 100 


47 ±2 


42 


204.28258 


-29.887080 


11 ±4 


19 ±4 


314 ±52 


720 ± 220 


340 ± 240 


- 


22 ±3 


43 


204.21262 


-29.876963 


34 ±4 


32 ±5 


884 ± 92 


1350 ±310 


660 ±310 


170 ± 160 


23 ±3 


44 


204.27939 


-29.857758 


10 ±4 


20 ±4 


312 ±60 


1360 ± 540 


630 ± 320 


390 ± 200 


25 ±3 


45 


204.28841 


-29.857795 


170 ± 10 


120 ±10 


3080 ± 210 


4000 ± 650 


1220 ± 470 


300 ± 200 


72 ±6 


46 


204.27239 


-29.821631 


17 ±4 


17 ±5 


424 ± 62 


1140 ±270 


770 ± 290 


370 ±180 


32 ±2 


47 


204.23190 


-29.881713 


27 ±5 


12 ±5 


724 ± 84 


2010 ± 680 


950 ± 350 


690 ± 220 


21 ±2 


48 


204.17584 


-29.875390 


37 ±4 


12 ±2 


266 ± 38 


1530 ± 230 


1180 ±270 


480 ± 110 


8±4 


49 


204.26047 


-29.832599 


25 ±3 


20 ±2 


374 ± 42 


1350 ± 280 


750 ± 380 


350 ± 190 


10 ±2 


50 


204.26316 


-29.827777 


92 ±9 


37 ±4 


674 ± 82 


1740 ± 370 


1410 ± 490 


210 ± 170 


34 ±2 


51 


204.26779 


-29.925076 


28 ±3 


19 ±2 


428 ± 34 


1550 ± 190 


810 ± 190 


400 ± 120 


38 ±3 


52 


204.28075 


-29.909675 


47 ±4 


37 ±3 


605 ± 44 


1400 ± 220 


610 ± 180 


250 ± 140 


77 ±3 


53 


204.32062 


-29.886212 


14 ±2 


12 ± 1 


223 ± 19 


410 ± 110 


450 ± 170 


270 ± 110 


23.5 ± 0.8 


54 


204.23495 


-29.885641 


15 ±5 


18 ±6 


294 ± 64 


780 ± 580 


280 ± 260 


- 


- 


55 


204.28509 


-29.878952 


26 ±5 


17 ±4 


539 ± 63 


1010 ±250 


600 ± 310 


240 ± 190 


22 ±4 


56 


204.27136 


-29.804832 


41 ±3 


27 ±2 


695 ± 43 


1450 ± 290 


590 ± 180 


250 ± 1 10 


39 ±2 


57 


204.24438 


-29.801254 


8±2 


10 ±2 


213 ±31 


280 ± 140 


360 ± 150 


200 ± 1 10 


26 ±2 


58 


204.22839 


-29.926435 


22 ±2 


9± 1 


287 ± 20 


1120 ± 170 


700 ± 190 


370 ± 110 


10 ± 0.7 


59 


204.17916 


-29.878426 


25 ±3 


9±2 


311 ±42 


1380 ± 240 


890 ± 240 


500 ± 110 


36 ±3 


60 


204.21423 


-29.871364 


20 ±4 


15 ±4 


296 ± 66 


800 ± 330 


860 ± 340 


- 


- 


61 


204.26419 


-29.850760 


17 ±4 


20 ±7 


515 ±90 


1050 ± 570 


710 ± 340 


220 ± 170 


51 ±3 


62 


204.26217 


-29.847377 


20 ±5 


18 ±6 


432 ± 88 


1230 ± 560 


640 ± 340 


190 ± 150 


86 ±6 


63 


204.30113 


-29.838620 


15 ±3 


9± 1 


265 ± 40 


800 ± 280 


640 ± 240 


350 ± 170 


14 ±2 
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Table D 1 - continued from previous page 



RA 
[deg] 



8 
[deg] 



[raJy] 



F24jjm 

[mJy] 



F70(im 

[mJy] 



ri60)im 
[mJy] 



F250m"i 
[mJy] 



[mJy] 



[/dy] 



64 


204.21281 


-29.835625 


65 


204.23465 


-29.825908 


66 


204.27451 


-29.895053 


67 


204.20733 


-29.869728 


68 


204.28385 


-29.911307 


69 


204.27882 


-29.887547 


70 


204.29993 


-29.850675 


71 


204.21230 


-29.881552 


72 


204.22276 


-29.853089 


73 


204.21339 


-29.841575 


74 


204.29638 


-29.826215 


75 


204.24362 


-29.804870 


76 


204.28500 


-29.866274 


77 


204.30424 


-29.860591 


78 


204.28876 


-29.852195 


79 


204.28907 


-29.846556 


80 


204.19144 


-29.880882 


81 


204.32204 


-29.864861 


82 


204.20639 


-29.859928 


83 


204.19682 


-29.818118 


84 


204.29686 


-29.926976 


85 


204.29044 


-29.908584 


86 


204.24101 


-29.839948 


87 


204.27353 


-29.917803 


88 


204.23122 


-29.914030 


89 


204.23179 


-29.802769 


90 


204.20338 


-29.873538 



24 ±3 


16 ±2 


365 ± 45 


890 ± 200 


540 ± 260 


310 ± 160 


46 ±2 


14 ±3 


6±4 


201 ± 50 


800 ± 240 


670 ±310 


250 ± 170 


17 ±2 


24 ±6 


11 ±3 


423 ± 62 


1180 ±260 


820 ± 260 


330 ± 180 


35 ±2 


17 ±4 


13 ±4 


387 ± 68 


450 ± 240 


570 ± 240 


200 ± 150 


26 ±4 


9±2 


6± 1 


194 ± 20 


390 ± 110 


290 ± 140 


150 ± 120 


5±2 


23 ±6 


11 ±3 


290 ± 49 


1590 ± 270 


1480 ± 360 


690 ± 220 


26 ±2 


10 ±3 


8±2 


195 ± 39 


410 ± 270 


330 ± 270 


180 ± 150 


25 ±2 


19 ±4 


7±5 


285 ± 70 


1200 ± 280 


740 ± 300 


300 ± 190 


4±3 


21 ±4 


24 ±7 


407 ± 74 


1360 ± 380 


750 ± 380 


440 ± 190 


12 ±3 


18 ±3 


26 ±4 


308 ± 46 


530 ± 220 


700 ± 320 


170 ± 140 


30 ±2 


11 ±4 


7±2 


157 ± 40 


570 ± 270 


480 ± 180 


300 ± 140 


9±4 


14 ±3 


9±2 


215 ±30 


680 ± 160 


530 ± 200 


350 ± 140 


23 ±2 


12 ±4 


16 ±3 


271 ±49 


1180 ±330 


680 ± 350 


300 ± 190 


19 ±2 


12 ±2 


5± 1 


185 ±28 


750 ± 200 


600 ± 260 


220 ± 150 


14 ± 1 


13 ±3 


4±3 


286 ± 49 


890 ± 280 


320 ± 300 


- 


16 ±3 


12 ±2 


6±3 


161 ± 33 


480 ± 260 


560 ± 350 


- 


20 ±3 


10 ±2 


5± 1 


135 ±20 


530 ± 190 


480 ± 200 


330 ± 120 


3±2 


10 ± 1 


4.6 ± 0.8 


110 ± 12 


410 ± 150 


230 ± 150 


110 ±72 


9±1 


14 ±4 


13±3 


224 ± 47 


570 ± 270 


350 ± 290 


430 ± 170 


16 ±3 


15 ±1 


8.5 ± 0.7 


146 ± 12 


584 ± 92 


340 ± 120 


180 ± 73 


12.8 ±0.5 


16 ±2 


13 ±1 


95 ±11 


400 ± 120 


290 ± 130 


145 ± 81 


21 ± 1 


10 ±2 


7 ±0.8 


94 ± 14 


309 ± 93 


310 ± 140 


220 ± 130 


3±1 


10 ±4 


5±2 


113 ±35 


840 ± 260 


660 ± 370 


340 ± 180 


34 ±2 


18 ±2 


12 ±2 


207 ± 30 


840 ± 170 


570 ± 180 


260 ± 120 


60 ±3 


8± 1 


3.4 ± 0.7 


96 ±14 


540 ± 170 


220 ± 160 


- 


5.5 ±0.7 


13 ±3 


3± 1 


134 ±29 


550 ± 180 


390 ± 140 


400 ± 120 


2±2 


11 ±4 


30 ±4 


148 ± 57 


1050 ± 260 


530 ± 250 


190 ± 150 


9±3 



Table D2: HB measurement SED fitting parameters 





XUV 


Xcol 


Tdust 
[K] 


Mdust 

[10 4 M e ] 


F 2 4 


Ldust 

[10 4() erg/s] 


Xfit 


1 


1.3±0.7 


9±2 


22.9±0.8 


16±5 


0.8±0.2 


47 ±1 


0.6 


2 


0.3±0.1 


9±1 


19.6±0.1 


18±3 


0.9±0.1 


28.7±0.6 


5.0 


3 


1.4±0.6 


6±2 


22.0±0.6 


20±5 


0.9±0.2 


45 ±2 


1.2 


4 


2±1 


8±3 


24±1 


12±5 


0.7±0.2 


35±2 


2.0 


5 


2±1 


4.±1 


23.6±0.7 


14±3 


0.4±0.2 


32 ±2 


1.0 


6 


3±2 


6±2 


26±1 


7±3 


0.6±0.1 


31±1 


0.6 


7 


7±4 


8±3 


31±3 


2±1 


0.6±0.2 


24±1 


1.1 


8 


1.3±0.7 


9±1 


23.1±0.6 


11±4 


0.8±0.2 


25 ±1 


6.6 


9 


2±1 


6±3 


22.9±0.7 


30±10 


0.6±0.1 


52 ±2 


4.5 


10 


2±1 


2±1 


21±1 


22±9 


0.7±0.2 


35±1 


2.5 


11 


1.4±0.6 


6±2 


21.8±0.4 


14±3 


0.7±0.2 


23 ±1 


0.5 


12 


1.4±0.6 


8±2 


22.6±0.4 


22±6 


0.6±0.1 


40.6 ±0.9 


0.4 


13 


3±2 


3±2 


23.2±0.3 


12±3 


0.4±0.1 


27.1 ±0.8 


0.4 


14 


0.6±0.2 


9±2 


21.6±0.6 


27±8 


0.5±0.2 


31 ±2 


1.9 


15 


0.4±0.3 


6±4 


18.2±0.8 


30±10 


0.7±0.2 


14.2 ±0.8 


2.3 


16 


1.4±0.6 


6±2 


21.6±0.2 


39±7 


0.8±0.1 


74.1±0.5 


2.9 


17 


2±1 


3.±1 


22.9±0.7 


4.4±0.8 


0.6±0.1 


10.3±0.3 


0.9 


18 


0.3±0.2 


9±2 


20±1 


13±6 


0.9±0.1 


14.6 ±0.8 


1.5 


19 


1.2±0.8 


7±3 


22±1 


8±4 


0.7±0.3 


12±1 


5.8 


20 


1.3±0.8 


7±3 


21.6±0.7 


13±4 


0.6±0.2 


16.2 ±0.8 


0.2 


21 


0.5±0.2 


9±2 


20.8±0.5 


17±4 


0.6±0.2 


17 ±1 


1.7 


22 


2±1 


4±2 


23.3±0.3 


8±2 


0.4±0.2 


18 ±1 


0.9 


23 


2±2 


6±4 


23±2 


5±3 


0.4±0.4 


7±1 


1.8 


24 


1.3±0.7 


3±2 


19.9±0.4 


20±5 


0.13±0.1 


12.9 ±0.8 


0.9 


25 


3±2 


4±2 


23.7±0.7 


14±4 


0.5±0.1 


37 ±1 


1.3 


26 


l.Oil.O 


7±4 


20±2 


9±6 


0.8±0.2 


10±1 


5.4 


27 


6±4 


5±3 


29±3 


1.5±1.0 


0.5±0.2 


8.8± 0.8 


0.6 


28 


1.4±0.6 


4±2 


20.8±0.4 


12±3 


0.13±0.1 


10.6 ±0.5 


0.7 


29 


0.5±0.3 


8±3 


20.2±0.8 


9±4 


0.7±0.2 


7.8 ±0.6 


1.1 
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Xuv 


,Ycol 


Tdust 
[K] 


M dust 
[10 4 M ] 


F24 


Ldust 

[10 4() erg/s] 


Xfit 


30 


0.6±0.4 


7±4 


20±1 


6±3 


0.9±0.2 


6.6± 0.5 


1.4 


31 


3±2 


1.2±0.9 


22±1 


9±4 


0.4±0.2 


17±1 


1.4 


32 


3±3 


7±4 


26±2 


2±1 


0.5±0.3 


6.8±0.7 


0.6 


33 


1.2±0.8 


7±3 


21.5±0.8 


5±2 


0.6±0.2 


5.9 ±0.4 


2.3 


34 


0.6±0.4 


8±3 


21±1 


8±4 


0.6±0.2 


7.3± 0.7 


0.5 


35 


0.2±0.1 


8±2 


18=1=1 


20±10 


0.6±0.3 


7.5 ±0.6 


4.5 


36 


3±2 


3±2 


22.8±0.7 


9±3 


0.4±0.1 


17.9± 0.6 


1.3 


37 


0.3=1=0.1 


9±1 


19.6±0.1 


14±3 


0.5±0.3 


9.1 ±0.5 


4.0 


38 


1.2±0.8 


7±3 


21.6±0.7 


16±6 


0.6±0.2 


21 ±1 


0.9 


39 


6±5 


3±3 


25±3 


2±2 


0.5±0.4 


7±1 


0.8 


40 


2±1 


4±2 


22±1 


6±3 


0.3±0.2 


7.2±0.5 


0.6 


41 


2±1 


3±2 


21.3=1=0.8 


6±3 


0.3±0.1 


7.2± 0.6 


0.9 


42 


1.2±0.9 


7±3 


22±2 


4±3 


0.7±0.3 


6.0± 0.7 


1.4 


43 


2±1 


6±2 


24.7±0.5 


5±1 


0.4±0.2 


14.1 ±0.6 


0.4 


44 


0.3±0.2 


9±2 


18.9=1=0.8 


13±6 


0.8±0.3 


7.5 ±0.8 


1.6 


45 


8±3 


2±1 


27=1=1 


8±3 


0.16±0.1 


48±2 


2.9 


46 


0.5±0.3 


8±2 


20.7±0.9 


12±5 


0.5±0.4 


9.0 ±0.8 


0.4 


47 


0.6±0.3 


9±1 


21.5±0.9 


15±6 


0.11±0.1 


13.6±0.9 


2.1 


48 


0.5±0.2 


4±2 


18.2=1=0.3 


25±6 


0.04±0.03 


9.8 ±0.3 


3.9 


49 


0.6±0.3 


5±2 


19.9±0.6 


14±5 


0.41±0.1 


9.6 ±0.5 


0.1 


50 


4±2 


0.5±0.4 


22.1 ±0.7 


10±3 


0.05±0.04 


17 ±1 


4.4 


51 


0.55±0.1 


6±2 


20.03±0.06 


15±2 


0.3±0.2 


10.7±0.3 


0.6 


52 


2±1 


2±1 


22.1 ±0.7 


7±1 


0.43±0.1 


13.3 ±0.5 


1.1 


53 


2±1 


3±2 


21=1=1 


4±2 


0.5±0.2 


4.8±0.4 


4.3 


54 


4±3 


6±5 


23±3 


5 ±4 


0.6±0.4 


6±1 


0.2 


55 


2±1 


7±4 


23=1=1 


6±3 


0.3±0.3 


10±1 


0.6 


56 


2±1 


3.±1 


22.9±0.7 


6.9±0.9 


0.21±0.08 


12.9±0.2 


1.0 


57 


1.2±0.9 


7±3 


22±2 


3±2 


0.6±0.3 


3.9± 0.5 


2.9 


58 


0.5±0.2 


7±2 


19.8±0.2 


12±2 


0.03±0.02 


7.45±0.04 


3.1 


59 


0.4±0.2 


7±3 


19.2±0.4 


19±6 


0.06±0.05 


9.0 ±0.5 


2.0 


60 


1.2±0.9 


5±3 


20±1 


9±4 


0.4±0.4 


7.1± 1.0 


1.5 


61 


1.2±0.9 


8±2 


22±2 


8±5 


0.5±0.4 


9±1 


0.3 


62 


1.2±0.9 


7±4 


22±1 


8±4 


0.5±0.4 


9±1 


0.4 


63 


0.5±0.3 


8±3 


20.1 ±0.9 


10±4 


0.2±0.2 


6.4 ±0.6 


0.5 


64 


1.3±0.7 


5±2 


21.1=1=0.6 


8±2 


0.3±0.2 


8.1 ±0.6 


0.6 


65 


0.5±0.4 


7±3 


19=1=1 


12±7 


0.3±0.3 


5.5±0.8 


0.5 


66 


0.6±0.4 


7±3 


21.1=1=0.8 


11±4 


0.14±0.1 


9.1 ±0.8 


0.8 


67 


3±2 


6±4 


24±2 


4±2 


0.4±0.4 


6.5±0.9 


2.2 


68 


1.2±0.8 


7±3 


21.8=1=0.8 


3±1 


0.3±0.3 


3.7 ±0.4 


0.7 


69 


0.20±0.1 


8±2 


18.1=1=0.4 


33±8 


0.12±0.1 


10.2±0.7 


1.0 


70 


2±1 


6±4 


22±2 


5±4 


0.4±0.4 


4.1 ±0.7 


0.3 


71 


0.5±0.3 


7±3 


19.7±0.9 


13±6 


0.2±0.2 


7.6 ±0.8 


0.5 


72 


0.4±0.3 


7±3 


19=1=1 


17±8 


0.6±0.4 


10± 1 


0.1 


73 


2±1 


3±3 


20±2 


8±5 


0.8±0.2 


7.0±0.6 


1.5 


74 


0.4±0.3 


7±4 


19=1=1 


11=1=6 


0.4±0.4 


4.6± 0.7 


0.5 


75 


0.6±0.4 


7±4 


19.7±0.9 


9±4 


0.3±0.3 


5.5 ±0.5 


1.3 


76 


0.3±0.2 


8±2 


18.9=1=0.8 


15±6 


0.6±0.3 


7.2 ±0.7 


0.7 


77 


0.4±0.2 


8±3 


19.6±0.8 


9±4 


0.09±0.08 


4.9 ±0.5 


0.8 


78 


1.2±0.9 


8±2 


22±1 


6±3 


0.12±0.1 


5.3±0.6 


2.0 


79 


1.1 ±0.9 


6±4 


20±1 


7±4 


0.3±0.3 


4.0± 0.7 


0.6 


80 


0.3±0.2 


7±3 


18.8±0.8 


10±4 


0.2±0.2 


4.1 ±0.5 


1.6 


81 


1.3±0.8 


4±2 


20.1 ±0.5 


4±2 


0.10±0.09 


3.0 ±0.2 


0.2 


82 


l.Oil.O 


6±4 


19±2 


10±7 


0.5±0.3 


5.6± 0.8 


2.0 


83 


1.3=1=0.8 


3±2 


19.5±0.4 


7±2 


0.21±0.1 


4.3 ±0.2 


0.5 


84 


0.6±0.4 


0.4±0.3 


16±1 


12±6 


0.6±0.2 


3.7 ±0.2 


0.9 


85 


1.2±0.9 


3±2 


19=1=1 


6±3 


0.4±0.2 


2.8 ±0.3 


1.3 


86 


0.3±0.2 


7±3 


17±1 


20±10 


0.2±0.2 


4.5 ±0.8 


0.5 


87 


0.6±0.3 


4±2 


19.2±0.7 


11±4 


0.3±0.2 


6.1 ±0.4 


0.1 


88 


1.2±0.9 


5±3 


19.9±1.0 


4±2 


0.07±0.06 


2.5 ±0.3 


2.8 


89 


0.4±0.3 


7±4 


18.6±0.9 


11=1=5 


0.06±0.05 


4.0 ±0.4 


6.3 


90 


0.20±0.1 


6±4 


16±1 


20±10 


0.8±0.2 


6.6 ±0.7 


6.6 
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Table D3: LB measurement source fluxes and positions 





RA 


5 


rSfim 


P24,fj.m 


F70^m 


Fl60(im 


r25Qfi.m 


r350^m 


F(H a ) 




[deg] 


[deg] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[MJy] 


1 


204.27867 


-29.826725 


132 ± 9.0 


302 ± 17 


3.17 ±0.17 


5.43 ± 0.55 


2.63 ± 0.67 


1.37 ±0.34 


70 ±7.1 


2 


204.26433 


-29.900421 


63 ± 6.9 


197 ± 12 


1.69 ±0.11 


2.59 ± 0.40 


1.36 ±0.54 


0.49 ± 0.32 


52 ± 2.2 


3 


204.22128 


-29.858999 


163 ± 12 


283 ± 18 


3.18 ±0.20 


5.03 ± 0.84 


2.25 ± 0.72 


0.70 ± 0.33 


133 ± 3.2 


4 


204.28271 


-29.854724 


114 ± 10. 


166 ± 11 


2.48 ±0.16 


4.40 ± 0.67 


1.93 ±0.72 


0.87 ± 0.35 


83 ±3.7 


5 


204.23187 


-29.830846 


142 ± 10. 


127 ± 8.7 


2.41 ±0.15 


4.04 ± 0.49 


1.70 ±0.57 


0.65 ± 0.28 


106 ± 3.9 


6 


204.22652 


-29.884745 


126 ±9.8 


192 ± 14 


2.49 ±0.17 


2.93 ±0.71 


1.26 ±0.58 


0.48 ±0.31 


205 ± 3.6 


7 


204.22894 


-29.885796 


84 ± 8.3 


144 ± 12 


2.32 ±0.17 


1.87 ±0.68 


0.87 ± 0.59 


- 


131 ±2.9 


8 


204.22172 


-29.882527 


65 ± 7.6 


146 ± 11 


2.08 ±0.15 


3.38 ± 0.65 


2.64 ± 0.80 


1.34 ±0.37 


79 ±3.1 


9 


204.21930 


-29.863860 


268 ± 17 


386 ± 25 


4.94 ±0.31 


7.7 ± 1.2 


2.60 ± 0.83 


1.57 ±0.44 


238 ± 3.7 


10 


204.22046 


-29.880233 


175 ± 14 


264 ± 19 


2.29 ±0.15 


5.24 ± 0.80 


3.18 ±0.85 


1.34 ±0.37 


153 ± 3.7 


11 


204.21971 


-29.852660 


93 ± 8.5 


106 ± 11 


1.56 ±0.13 


3.21 ±0.51 


1.54 ±0.69 


0.42 ±0.31 


91 ±3.6 


12 


204.26957 


-29.849437 


159 ± 13 


190 ± 13 


3.09 ± 0.21 


5.41 ±0.98 


2.16 ±0.87 


1.05 ±0.40 


210 ±3.7 


13 


204.28487 


-29.869891 


144 ± 13 


125 ±9.6 


1.79 ±0.13 


3.34 ± 0.54 


1.37 ±0.67 


0.46 ±0.31 


101 ±3.1 


14 


204.22800 


-29.882334 


116 ±9.4 


132 ± 11 


2.17 ±0.16 


4.14 ±0.86 


3.37 ±0.80 


1.27 ±0.40 


34 ± 2.5 


15 


204.28634 


-29.858535 


70 ± 8.4 


106 ± 11 


0.796 ± 0.080 


2.93 ± 0.52 


1.26 ±0.67 


0.69 ± 0.33 


20 ± 3.6 


16 


204.28083 


-29.852275 


260 ± 20 


422 ± 32 


5.04 ± 0.34 


8.1 ± 1.1 


5.5 ± 1.4 


1.92 ±0.46 


193 ± 4.5 


17 


204.16814 


-29.855863 


47 ± 2.9 


48 ± 3.0 


0.665 ± 0.037 


1.39 ±0.19 


0.85 ± 0.22 


0.45 ±0.12 


62 ± 0.93 


18 


204.22554 


-29.844806 


51 ±5.8 


93 ± 7.2 


0.918 ± 0.085 


2.19 ±0.45 


0.73 ± 0.62 


- 


83 ±3.2 


19 


204.21463 


-29.883342 


58 ± 7.4 


91 ± 9.7 


1.06 ±0.11 


1.70 ±0.42 


2.07 ± 0.73 


0.82 ± 0.32 


55 ±3.6 


20 


204.28397 


-29.881946 


91 ± 10.0 


76 ± 7.8 


1.200 ±0.10 


2.47 ± 0.44 


1.11 ±0.60 


0.31 ±0.29 


76 ± 3.7 


21 


204.20849 


-29.878634 


77 ±8.1 


79 ± 7.3 


1.21 ±0.12 


3.75 ± 0.52 


2.03 ± 0.74 


1.14 ±0.36 


47 ±4.1 


22 


204.26997 


-29.823576 


97 ±9.1 


81 ±8.0 


1.39 ±0.12 


3.10 ±0.45 


1.94 ±0.62 


0.78 ± 0.28 


135 ±7.3 


23 


204.29186 


-29.858028 


42 ± 5.9 


38 ± 7.0 


0.619 ± 0.076 


1.00 ±0.40 


0.93 ± 0.65 


0.51 ±0.33 


98 ±4.7 


24 


204.25244 


-29.905010 


66 ± 9.6 


40 ± 4.0 


0.606 ± 0.072 


1.62 ±0.35 


0.86 ± 0.48 


0.29 ± 0.29 


54 ± 2.3 


25 


204.21223 


-29.844608 


186 ± 12 


165 ± 11 


2.82 ±0.18 


5.53 ±0.58 


2.50 ±0.61 


1.37 ±0.33 


184 ±4.1 


26 


204.21877 


-29.863393 


74 ± 7.2 


147 ± 13 


1.36 ±0.13 


2.17 ±0.74 


1.96 ±0.72 


- 


155 ±2.8 


27 


204.28831 


-29.849784 


50 ±5.9 


51 ±7.4 


0.849 ± 0.085 


0.97 ± 0.40 


0.80 ± 0.60 


0.65 ± 0.30 


67 ±4.2 


28 


204.24629 


-29.907226 


38 ± 6.5 


26 ± 3.3 


0.474 ± 0.071 


0.93 ± 0.32 


0.47 ± 0.45 


- 


49 ± 2.3 


29 


204.18835 


-29.879014 


32 ± 3.9 


32 ± 2.6 


0.473 ± 0.045 


1.21 ±0.30 


0.83 ± 0.46 


0.23 ±0.19 


104 ± 3.5 


30 


204.18091 


-29.872982 


22 ± 3.8 


41 ±3.6 


0.459 ± 0.049 


0.75 ± 0.28 


0.72 ± 0.42 


- 


71 ±4.5 


31 


204.23124 


-29.845585 


117 ± 11 


105 ± 8.6 


1.066 ±0.092 


2.37 ± 0.56 


0.78 ± 0.59 


0.31 ±0.30 


29 ± 2.5 


32 


204.29158 


-29.819528 


31 ±4.6 


32 ±3.1 


0.620 ± 0.070 


1.25 ±0.31 


0.93 ± 0.25 


0.81 ±0.21 


109 ± 4.3 


33 


204.22354 


-29.813422 


20 ± 3.3 


22 ± 2.2 


0.328 ± 0.036 


0.63 ± 0.23 


0.48 ± 0.39 


0.28 ± 0.23 


19 ±2.8 


34 


204.29732 


-29.831076 


49 ± 5.3 


40 ± 3.3 


0.709 ± 0.063 


2.64 ± 0.39 


2.04 ± 0.42 


1.12 ±0.27 


76 ± 3.3 


35 


204.27795 


-29.823031 


36 ±4.8 


51 ±6.0 


0.518 ±0.063 


3.38 ± 0.45 


1.79 ±0.57 


0.95 ± 0.30 


18 ±7.1 


36 


204.24819 


-29.810471 


84 ± 9.8 


67 ± 5.0 


1.110 ±0.086 


2.30 ± 0.37 


1.14 ±0.50 


0.42 ± 0.21 


87 ±7.0 


37 


204.26738 


-29.899060 


38 ±5.8 


46 ± 5.7 


0.586 ± 0.069 


1.80 ±0.40 


0.91 ±0.50 


- 


84 ± 2.0 


38 


204.23665 


-29.880081 


102 ± 12 


106 ± 8.4 


1.69 ±0.16 


4.3 ±1.0 


2.56 ± 0.70 


1.03 ±0.38 


36 ± 2.4 


39 


204.21163 


-29.866269 


42 ± 6.8 


39 ± 6.3 


0.616 ± 0.096 


1.59 ±0.38 


1.42 ±0.54 


0.72 ± 0.29 


42 ± 4.3 


40 


204.31140 


-29.843049 


46 ± 4.3 


29 ± 2.3 


0.458 ± 0.041 


1.92 ±0.32 


1.15 ±0.29 


0.68 ±0.19 


42 ± 15 


41 


204.25832 


-29.925232 


42 ± 3.7 


30 ± 2.2 


0.447 ± 0.032 


1.27 ±0.23 


0.72 ± 0.25 


0.39 ±0.16 


66 ±1.5 


42 


204.28258 


-29.887080 


35 ± 5.3 


41 ± 4.7 


0.561 ± 0.065 


1.95 ±0.32 


1.11 ±0.39 


0.82 ± 0.32 


46 ± 3.3 


43 


204.21262 


-29.876963 


86 ± 8.3 


93 ± 7.8 


1.64 ±0.14 


2.97 ± 0.50 


1.72 ±0.72 


0.77 ± 0.29 


65 ±3.8 


44 


204.27939 


-29.857758 


29 ± 5.6 


44 ± 6.3 


0.483 ± 0.076 


3.18 ±0.72 


1.25 ±0.62 


0.81 ±0.35 


35 ±3.5 


45 


204.28841 


-29.857795 


321 ± 28 


271 ± 26 


4.76 ± 0.34 


6.6 ±1.1 


1.59 ±0.88 


0.63 ± 0.34 


161 ±6.5 


46 


204.27239 


-29.821631 


57 ± 6.5 


51 ±6.7 


0.782 ± 0.085 


2.94 ± 0.44 


1.79 ±0.59 


0.90 ± 0.29 


130 ±7.0 


47 


204.23190 


-29.881713 


69 ± 7.6 


64 ±8.1 


1.22 ±0.12 


2.92 ± 0.83 


1.94 ±0.67 


1.09 ±0.38 


68 ± 2.2 


48 


204.17584 


-29.875390 


77 ± 7.8 


27 ± 3.2 


0.490 ± 0.053 


1.94 ±0.33 


1.24 ±0.50 


0.20 ±0.17 


22 ± 4.4 


49 


204.26047 


-29.832599 


46 ±5.1 


35 ± 2.9 


0.446 ± 0.054 


2.36 ± 0.40 


1.75 ±0.62 


0.80 ±0.31 


- 


50 


204.26316 


-29.827777 


140 ± 15 


66 ± 6.4 


0.87 ±0.11 


2.59 ±0.50 


2.08 ± 0.73 


0.42 ± 0.27 


43 ±3.1 


51 


204.26779 


-29.925076 


53 ±4.1 


38 ± 2.9 


0.635 ± 0.045 


2.16 ±0.26 


1.13 ±0.29 


0.49 ±0.17 


63 ±3.3 


52 


204.28075 


-29.909675 


80 ± 6.2 


64 ±4.1 


0.854 ± 0.057 


2.38 ±0.31 


1.14 ±0.30 


0.52 ± 0.21 


135 ± 3.2 


53 


204.32062 


-29.886212 


26 ± 2.9 


19 ± 1.9 


0.294 ± 0.029 


0.72 ±0.17 


0.84 ±0.31 


0.35 ±0.17 


36± 1.1 


54 


204.23495 


-29.885641 


38 ±6.8 


42 ±7.1 


0.551 ±0.088 


1.33 ±0.70 


0.88 ±0.51 


0.32 ± 0.27 


15 ±2.6 


55 


204.28509 


-29.878952 


65 ± 8.0 


41 ± 6.7 


0.838 ± 0.089 


1.55 ±0.41 


0.70 ± 0.63 


- 


61 ±3.9 


56 


204.27136 


-29.804832 


51 ±9.2 


36 ± 3.6 


0.753 ± 0.097 


1.31 ±0.41 


0.60 ± 0.41 


0.25 ±0.16 


46 ± 6.9 


57 


204.24438 


-29.801254 


15 ±3.8 


17 ±2.5 


0.257 ± 0.043 


0.73 ± 0.26 


0.86 ± 0.48 


0.36 ±0.19 


18 ±7.1 


58 


204.22839 


-29.926435 


37 ± 3.0 


18 ± 1.5 


0.385 ± 0.030 


1.44 ±0.22 


0.81 ±0.28 


0.30 ±0.16 


15 ± 0.82 


59 


204.17916 


-29.878426 


54 ± 5.6 


27 ± 3.4 


0.561 ± 0.060 


1.79 ±0.35 


1.04 ±0.48 


0.25 ±0.19 


139 ± 3.6 


60 


204.21423 


-29.871364 


40 ±6.1 


31 ±5.4 


0.440 ± 0.082 


1.58 ±0.48 


1.52 ±0.68 


0.43 ± 0.29 


- 


61 


204.26419 


-29.850760 


53 ± 6.9 


66 ± 9.3 


1.01 ±0.12 


2.89 ± 0.82 


1.13 ±0.76 


0.71 ± 0.37 


98 ±4.0 


62 


204.26217 


-29.847377 


53 ± 7.6 


48 ± 7.5 


0.79 ±0.12 


1.75 ±0.73 


0.80 ± 0.78 


0.37 ± 0.30 


175 ±6.1 


63 


204.30113 


-29.838620 


47 ± 4.4 


28 ± 2.3 


0.514 ± 0.052 


2.46 ± 0.39 


2.09 ± 0.44 


1.18 ±0.28 


39 ± 2.0 
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Table D3 - continued from previous page 



RA 
[deg] 



<5 
[deg] 



[mJy] 



F24jjm 

[mJy] 



F70(im 

[mJy] 



[mJy] 



[mJy] 



[mJy] 



HH a ) 
[fjJy] 



64 


204.21281 


-29.835625 


65 


204.23465 


-29.825908 


66 


204.27451 


-29.895053 


67 


204.20733 


-29.869728 


68 


204.28385 


-29.911307 


69 


204.27882 


-29.887547 


70 


204.29993 


-29.850675 


71 


204.21230 


-29.881552 


72 


204.22276 


-29.853089 


73 


204.21339 


-29.841575 


74 


204.29638 


-29.826215 


75 


204.24362 


-29.804870 


76 


204.28500 


-29.866274 


77 


204.30424 


-29.860591 


78 


204.28876 


-29.852195 


79 


204.28907 


-29.846556 


80 


204.19144 


-29.880882 


81 


204.32204 


-29.864861 


82 


204.20639 


-29.859928 


83 


204.19682 


-29.818118 


84 


204.29686 


-29.926976 


85 


204.29044 


-29.908584 


86 


204.24101 


-29.839948 


87 


204.27353 


-29.917803 


88 


204.23122 


-29.914030 


89 


204.23179 


-29.802769 


90 


204.20338 


-29.873538 



54 ±5.1 


33 ±3.1 


0.585 ± 0.058 


1.80 ±0.29 


1.07 ±0.40 


0.60 ± 0.23 


61 ±3.2 


35 ± 4.9 


20 ± 4.4 


0.369 ± 0.064 


1.31 ±0.33 


0.89 ± 0.48 


0.27 ± 0.24 


- 


46 ± 9.2 


22 ±4.1 


0.54 ±0.11 


1.32 ±0.42 


0.84 ± 0.54 


- 


49 ± 2.0 


40 ± 6.3 


32 ± 4.9 


0.651 ±0.087 


1.57 ±0.39 


1.36 ±0.50 


0.73 ± 0.29 


60 ± 4.2 


22 ± 2.7 


17 ± 1.7 


0.296 ± 0.026 


1.05 ±0.17 


0.84 ± 0.27 


0.43 ±0.19 


11 ± 1.8 


72 ± 10. 


28 ± 6.5 


0.441 ± 0.069 


2.16 ±0.44 


1.62 ±0.66 


0.61 ± 0.32 


116 ±2.9 


38 ±4.3 


20 ± 2.5 


0.399 ± 0.050 


1.80 ±0.36 


2.00 ± 0.50 


1.05 ±0.27 


44 ± 2.4 


64 ± 7.5 


48 ± 7.0 


0.699 ± 0.095 


2.29 ± 0.44 


1.74 ±0.71 


0.78 ± 0.32 


42 ± 3.5 


57 ± 6.5 


65 ±9.1 


0.750 ± 0.099 


2.44 ± 0.54 


1.06 ±0.63 


0.43 ±0.31 


15 ±3.2 


61 ±5.2 


63 ± 5.2 


0.685 ± 0.065 


2.01 ± 0.33 


1.63 ±0.51 


0.65 ± 0.23 


74 ± 3.4 


35 ±5.1 


22 ± 3.0 


0.346 ± 0.050 


1.78 ±0.35 


1.39 ±0.32 


1.18 ±0.27 


29 ± 3.7 


32 ± 5.3 


21 ±3.1 


0.346 ±0.051 


1.34 ±0.29 


1.16 ±0.52 


0.57 ± 0.22 


36 ± 7.0 


38 ±5.7 


45 ± 7.2 


0.549 ± 0.075 


2.46 ± 0.52 


0.87 ± 0.66 


0.34 ± 0.30 


52 ± 2.5 


32 ± 3.4 


15 ±2.0 


0.325 ± 0.038 


1.80 ±0.30 


1.59 ±0.45 


0.77 ± 0.25 


22 ± 1.4 


32 ± 5.0 


23 ± 5.8 


0.503 ± 0.070 


1.52 ±0.43 


0.74 ± 0.59 


0.63 ± 0.30 


68 ± 3.5 


29 ± 4.3 


19 ±5.8 


0.306 ± 0.053 


0.74 ± 0.39 


0.97 ± 0.66 


0.60 ± 0.30 


87 ± 3.7 


23 ± 3.6 


14 ± 1.9 


0.221 ±0.032 


0.94 ± 0.29 


0.77 ± 0.46 


0.32 ± 0.22 


12 ±3.3 


24 ±2.1 


12 ± 1.1 


0.200 ±0.018 


1.19 ±0.20 


0.97 ± 0.25 


0.66 ±0.14 


19 ± 1.2 


27 ± 6.2 


26 ± 4.2 


0.340 ± 0.064 


1.05 ±0.39 


1.20 ±0.64 


0.85 ± 0.29 


34 ±3.8 


26 ± 1.7 


16± 1.1 


0.225 ± 0.017 


1.00 ±0.13 


0.70 ±0.19 


0.43 ±0.12 


36 ± 0.54 


31 ±3.4 


24 ± 1.9 


0.170 ±0.017 


0.62 ±0.17 


0.44 ± 0.23 


0.18 ±0.13 


39 ± 1.4 


22 ± 3.2 


14 ± 1.6 


0.133 ±0.023 


0.71 ±0.16 


0.82 ± 0.29 


0.47 ± 0.20 


7.5 ± 1.3 


37 ± 6.2 


13 ±2.8 


0.239 ± 0.055 


1.45 ±0.49 


0.90 ± 0.85 


0.36 ± 0.34 


167 ± 2.6 


44 ± 4.0 


31 ±2.7 


0.429 ± 0.043 


1.58 ±0.25 


1.05 ±0.31 


0.44 ±0.18 


190 ± 2.6 


18 ±2.2 


7.6 ± 1.0 


0.178 ±0.025 


0.84 ± 0.22 


0.35 ± 0.25 


- 


7.30 ± 0.75 


21 ± 4.7 


7.1 ± 1.9 


0.123 ±0.051 


0.58 ± 0.28 


0.64 ± 0.39 


0.36 ± 0.21 


- 


30 ± 5.7 


74 ± 7.4 


0.380 ± 0.078 


1.74 ± 0.39 


1.15 ±0.53 


0.66 ± 0.30 


25 ± 3.7 



Table D4: LB measurement SED fitting parameters 





XUV 


Xcol 


Tdust 
[K] 


Mdusi 

[10 4 M ] 


F24 


Ldust 

[10 40 erg/s] 


Xfit 


1 


0.7±0.3 


4±2 


19.8±0.7 


40±20 


0.8±0.2 


67 ±2 


1.3 


2 


0.5±0.3 


4±3 


18.1±1.0 


30±10 


0.9±0.1 


36 ±1 


1.7 


3 


3±2 


2±1 


22.3±0.2 


20±5 


0.8±0.1 


65± 1 


1.9 


4 


2±1 


4±2 


22.1 ±0.9 


30±10 


0.7±0.1 


49 ±2 


1.4 


5 


3±1 


2.±1 


23.5±0.7 


17±2 


0.40±0.08 


44.8±0.5 


1.5 


6 


6±4 


2±1 


25±2 


8±5 


0.8±0.1 


45 ±3 


0.8 


7 


7±4 


4±2 


28±2 


4±2 


0.7±0.2 


36±2 


0.4 


8 


0.5±0.3 


9±1 


20.7±0.9 


30±10 


0.8±0.2 


41 ±2 


2.3 


9 


3±2 


2±1 


22.8±0.7 


30±10 


0.8±0.1 


97 ±2 


2.3 


10 


1.3±0.8 


1±1 


18±1 


70±30 


0.7±0.1 


59 ±2 


1.6 


11 


2±1 


2±1 


22.1 ±0.7 


15±2 


0.6±0.2 


32 ±1 


1.5 


12 


2±1 


3±1 


22.9±0.7 


23±4 


0.6±0.2 


58 ±3 


1.1 


13 


4±2 


1.2±0.8 


23±1 


14±5 


0.5±0.1 


38 ±2 


1.0 


14 


1.3±0.7 


5±2 


20.8±0.4 


37±9 


0.7±0.1 


46 ±1 


1.6 


15 


0.5±0.3 


1.1 ±0.9 


16.8±0.9 


50±20 


0.7±0.1 


24±1 


4.2 


16 


2±1 


3±2 


21±1 


50±20 


0.8±0.2 


106±5 


2.3 


17 


1.5±0.5 


3±2 


20.2±0.4 


14±3 


0.61±0.08 


15.6±0.2 


2.7 


18 


1.2±0.8 


3±2 


19±1 


20±8 


0.9±0.2 


21.6 ±0.9 


2.2 


19 


1.2±0.9 


3±2 


19±2 


30±10 


0.8±0.2 


23± 1 


2.8 


20 


3±2 


2±1 


22.8±0.7 


12±5 


0.4±0.2 


25 ±2 


1.5 


21 


0.6±0.3 


5±2 


19.5±0.6 


40±20 


0.5±0.1 


30 ±1 


0.4 


22 


2±1 


3±2 


21.3±0.8 


23±8 


0.5±0.2 


30± 1 


1.1 


23 


3±2 


2±2 


22±1 


8±5 


0.5±0.3 


12.7±0.9 


1.0 


24 


3±2 


1.2±0.9 


21.5±0.9 


11±4 


0.2±0.2 


15 ±1 


0.2 


25 


2±1 


2±1 


22.1 ±0.7 


32±4 


0.45±0.08 


58.1±0.7 


1.2 


26 


2±2 


1±1 


20±3 


20±10 


0.9±0.2 


30±2 


3.1 


27 


4±3 


3±2 


24±2 


6±3 


0.6±0.2 


16±1 


2.4 


28 


5 ±4 


2±1 


24±2 


4±3 


0.3±0.3 


10±1 


0.4 


29 


1.3±0.7 


4±2 


20.2±0.3 


11±3 


0.6±0.1 


11.2 ±0.6 


0.6 
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Table D4 - continued from previous page 





Xuv 


Xcol 


Tdust 
[K] 


Mdusi 

[10 4 M ] 


F24 


Ldust 

[10 40 erg/s] 


Xfit 


30 


2±1 


4±3 


20±2 


11±8 


0.9±0.2 


10.0±0.9 


0.6 


31 


3±2 


0.3±0.2 


21±1 


15±6 


0.5±0.2 


27± 1 


1.7 


32 


0.6±0.4 


7±3 


20.6±0.9 


15±6 


0.5±0.1 


13.4 ±0.8 


4.7 


33 


2±1 


3±2 


21±1 


6±4 


0.6±0.2 


7.2±0.7 


0.5 


34 


0.4±0.2 


6±2 


18.7±0.2 


40±10 


0.4±0.1 


20.1 ±0.9 


1.2 


35 


0.15±0.05 


9±2 


16.5±0.2 


80±20 


0.7±0.1 


18.9±0.9 


4.2 


36 


2±1 


2±1 


22.1 ±0.7 


13±2 


0.45±0.08 


23.7±0.3 


0.7 


37 


1.2±0.9 


4±3 


20±1 


19±9 


0.7±0.2 


15± 1 


1.2 


38 


1.3±0.8 


4±2 


20.6±0.6 


40±10 


0.57±0.1 


38 ±1 


0.3 


39 


1.2±0.9 


5±3 


20±1 


20±10 


0.5±0.2 


15± 1 


1.2 


40 


0.55±0.1 


3.±1 


18.8±0.3 


27±5 


0.3±0.1 


14.1±0.2 


0.3 


41 


2±1 


2±1 


20.3±0.8 


12±5 


0.4±0.1 


11.7 ±0.9 


1.3 


42 


0.4±0.3 


6±3 


18.7±0.6 


30±10 


0.6±0.2 


15.2 ±0.9 


0.6 


43 


2±1 


4±2 


22.1 ±0.9 


20±8 


0.5±0.1 


32±2 


1.1 


44 


0.20±0.1 


8±3 


17.1±0.9 


50±20 


0.7±0.2 


16 ±1 


2.8 


45 


8±2 


1.3±0.7 


26.4±0.8 


14±4 


0.4±0.2 


83 ±4 


2.7 


46 


0.5±0.3 


5±2 


19.0±0.6 


40±10 


0.5±0.2 


22 ±1 


0.3 


47 


1.2±0.8 


5±2 


20.7±0.7 


30±10 


0.5±0.2 


28 ±2 


0.6 


48 


2±1 


1.2±0.8 


20.8±0.3 


12±2 


0.02±0.02 


14.2±0.7 


7.5 


49 


0.4±0.2 


4±2 


17.8±0.1 


40±10 


0.4±0.1 


15.8±0.9 


0.3 


50 


3±2 


0.4±0.3 


21.3±1.0 


20±8 


0.09±0.08 


26 ±1 


2.0 


51 


1.4±0.6 


3±1 


20.2±0.3 


20±5 


0.33±0.1 


16.8±0.5 


1.1 


52 


2±1 


1.4±0.6 


21.0±0.1 


15±2 


0.40±0.08 


21.3±0.4 


1.7 


53 


1.4±0.6 


2±1 


19.7±0.2 


10±2 


0.4±0.2 


7.6 ±0.3 


2.4 


54 


2±1 


3±2 


20±2 


10±10 


0.6±0.2 


13±2 


0.2 


55 


6±4 


2±1 


25±1 


6±3 


0.2±0.2 


16±2 


1.0 


56 


3±2 


3±2 


23.8±0.9 


6±3 


0.3±0.2 


14±1. 


0.1 


57 


0.5±0.4 


7±4 


19±1 


14±8 


0.7±0.2 


6.6± 0.8 


1.0 


58 


1.4±0.6 


3.±1 


20.2±0.3 


13±3 


0.09±0.08 


10.5±0.3 


0.9 


59 


2±1 


3±2 


21.3±0.8 


13±6 


0.08±0.07 


14 ±1 


2.3 


60 


1.2±0.9 


3±2 


18.9±1.0 


20±10 


0.5±0.3 


13 ±1 


0.6 


61 


1.2±0.9 


6±3 


20±1 


20±10 


0.7±0.2 


23±2 


0.5 


62 


2±2 


3±2 


22±2 


12±7 


0.5±0.2 


17±2 


0.1 


63 


0.3±0.1 


5±1 


18.1±0.5 


47±6 


0.23±0.10 


17.1±0.3 


1.2 


64 


1.4±0.6 


3±2 


20.1 ±0.5 


20±5 


0.26±0.1 


15.7 ±0.8 


0.9 


65 


1.3±0.8 


3±2 


19.8±0.8 


15±6 


0.3±0.3 


10 ±1 


0.6 


66 


3±2 


3±2 


22.8±0.7 


6±2 


0.10±0.09 


11.3 ±0.9 


0.5 


67 


1.2±0.8 


6±3 


20.5±0.9 


17±7 


0.4±0.3 


15 ±1 


1.5 


68 


0.5±0.3 


5±2 


19.1±0.6 


16±6 


0.4±0.2 


8.2 ±0.6 


0.9 


69 


1.3±0.8 


2±1 


18.6±0.6 


30±10 


0.11 ±0.1 


16±1 


1.2 


70 


0.3±0.2 


6±3 


17.9±0.5 


40±20 


0.15±0.1 


13.7±0.9 


3.1 


71 


1.2±0.8 


3±2 


19.2±0.8 


30±10 


0.4±0.2 


19 ±1 


0.5 


72 


1.2±0.8 


3±2 


19.3±0.9 


30±10 


0.6±0.2 


20 ±1 


1.3 


73 


0.7±0.3 


2.±1 


18.6±0.7 


30±10 


0.6±0.2 


19.4±0.9 


0.7 


74 


0.3±0.2 


6±3 


17.6±0.5 


50±10 


0.3±0.2 


12.6 ±0.8 


2.7 


75 


0.6±0.4 


4±2 


19.0±0.9 


20±10 


0.3±0.2 


10.5 ±0.8 


0.7 


76 


0.5±0.3 


5±3 


18.6±1.0 


30±10 


0.7±0.2 


15 ±1 


2.7 


77 


0.3±0.2 


6±2 


18.1±0.4 


40±10 


0.07±0.06 


11.4±0.6 


1.5 


78 


1.2±0.9 


6±3 


20±1 


17±8 


0.4±0.4 


12± 1 


0.7 


79 


2±1 


3±2 


20±2 


10±10 


0.4±0.4 


8±1 


1.9 


80 


1.2±0.9 


3±2 


19±1 


14±7 


0.3±0.2 


7.1 ±0.8 


0.1 


81 


0.3±0.2 


5±2 


17.4±0.3 


30±10 


0.10±0.09 


8.0±0.5 


2.2 


82 


0.4±0.3 


6±4 


17±2 


40±30 


0.6±0.3 


10± 1 


2.1 


83 


0.7±0.3 


3±2 


18.8±0.4 


14±4 


0.27±0.10 


7.6 ±0.4 


1.2 


84 


0.6±0.4 


0.2±0.1 


17±1 


19±9 


0.5±0.2 


6.8±0.3 


2.6 


85 


0.4±0.3 


2±1 


16.5±0.8 


30±10 


0.4±0.2 


5.7 ±0.5 


1.1 


86 


1.2±0.9 


2±2 


18.7±0.8 


20±10 


0.07±0.06 


9±1 


1.3 


87 


1.3±0.7 


2±1 


19.4±0.5 


19±6 


0.36±0.1 


12.8±0.4 


0.3 


88 


1.2±0.8 


4±2 


19.8±0.8 


8±3 


0.07±0.06 


5.1 ±0.4 


1.8 


89 


1.1 ±0.9 


3±2 


18±2 


20±10 


0.10±0.09 


4.6 ±0.7 


1.4 


90 


0.2±0.1 


4±3 


14.9±0.7 


70±30 


0.9±0.2 


15 ±1 


7.4 



